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Abstract 


In the present thesis, the torque and the flux of a three-level induction motor are 
directly and independently controlled by using “Direct Torque and Flux Control (DTFC)” 
method. A three-level Synchronous Link Converter (SLC) is used as the front-end 
converter for drawing unity power factor sinusoidal current from utility. Hence the drive 
is utility friendly and capable of four quadrant operation. High performance induction 
motor drives require accurate independent torque and flux control over a wide range of 
speed range. DTFC is less complex and gives better control characteristics. The main 
difference between DTFC and traditional AC drive method controls is that with DTFC 
there is no separate voltage and frequency-controlled PWM modulator. 

The DTFC method gives attractive performance in terms of fast torque response, 
simple control scheme without coordinate transformation, and robustness against the 
motor parameter variation. It has been developed on the space vector approach, where the 
torque and flux of an induction motor can be controlled directly and independently. In the 
DTFC, the motor torque and flux are calculated from the primary variables and they are 
controlled directly and independently by selecting optimum inverter switch modes. The 
conventional DTFC suffers from flux demagnetization at low speed due to increased 
stator resistance drop. A modified control scheme has been developed in the present work 
to avoid this problem for better response at low speed condition. The real time algorithm 
of the drive system has been implemented with a PC, which uses PCL-208 data 
acquisition card for data transfer. 

Keywords: Direct Torque and Flux Control (DTFC), Synchronous Link Converter (SLC), 
Three-level inverter. Induction Motor, Look up table. Space vector modulation. 
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Chapter 1 


Introduction 


1.1 Background 

Motion control is required in large a number of industrial and domestic 
applications like transportation systems, rolling mills, paper machines, textile mills, 
machine tools, fans, pumps, robots, washing machines, etc. One major application of 
electric drives is in electric traction. Electric traction includes electric trains, electric 
buses, trams, battery driven and solar powered vehicles. 

In the field of electric traction, the ac electrification at commercial frequency 
began during the year 1955. The overhead catenary supply was fed at 25.0 kV, 50 Hz, 
single phase supply. Such a system needed power conversion on the Locomotives / Motor 
coaches for feeding the dc traction motors. Initially power conversion was done by means 
of mercury arc rectifiers. The maintenance cost of mercury arc rectifier was very high. In 
order to overcome this drawback, semiconductor diode rectifier was used to feed the dc 
traction motors. Such system has drawbacks such as, regeneration was not possible 
resulting in wastage of considerable wastage of power in the dynamic braking resistors, 
reduced life and increased maintenance of braking equipment. The ixiulti-stage phase 
controlled converters operating in sequence control and the pulse width modulated 



converters were used for obtaining the continuous control of tractive force/speed. 
Although these schemes ensured the better utilization of adhesion, the system has 
drawbacks due to the use of dc traction motors such as: the increased maintenance due to 
presence of commutator and brush gears, comparatively lower torque to weight ratio, 
lower operating speeds, lower efficiency compared with the ac traction motors. 

To overcome these drawbacks, the three phase ac motors have been employed. 
Such a drive has several important features like regenerative braking capability, very 
little maintenance, better utilization of available adhesion, higher average speed, higher 
torque-to-weight ratio, higher voltage and higher efficiency. Both synchronous and the 
asynchronous induction motors are used in the three phase ac motor traction drives. 

For almost a century, induction machine (IM) has been the work-horse of industry 
due to its robustness, low cost, and less maintenance. Before power electronics was 
introduced, induction motors were mainly used for essentially constant speed applications 
because of the unavailability of the variable-frequency voltage or current supply. The 
advancement of power electronics has made it possible to vary the frequency of the 
voltage or current supplies relatively easily, thus has extended the used of induction 
motor in variable speed drive applications; but due to the inherent coupling of flux and 
torque components in IM, it could not provide the torque performance comparable to the 
DC machine. 

The only effective way of producing a variable speed induction motor drive is to 
supply the induction motor with three phase voltages of variable frequency and variable 
amplitude. A variable frequency is required because the rotor speed depends on the speed 
of the rotating magnetic field provided by the stator. A variable voltage is required 
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because the motor impedance at low frequencies and consequently the current has to be 
limited by the means of reducing the supply voltage. 

It is in this area considerable research effort has been devoted with the aim to find 
even simpler methods of speed control for induction machines. One method, which has 
gained popularity over the last decade, is Direct Torque and Flux Control. 

Historically, several general controllers have been developed for induction motor control: 

■ Scalar controllers: Despite the fact that “Voltage-Frequency” (V/f) is the simplest 
controller, it is the most widespread, and is used in the majority of the industrial 
applications. It is known as scalar control [1] as it acts by imposing a constant relation 
between voltage and frequency without any phase angle control. The controller 
structure is very simple and it is normally used without speed feedback. In this 
control scheme, the performance of the machine improves in the steady state only, but 
the transient response is poor. However, this controller doesn’t achieve a good 
accuracy in both speed and torque responses, mainly due to the fact that the stator 
flux and torque are not directly controlled. 

■ Vector controllers: In 1971, Blaschke [2] proposed a scheme, which aims as the 
control of an induction motor like a separately excited dc motor, called field 
orientation control, vector control, or transvector control [3]. The idea of vector 
controlled AC motor drives is based on the excellent control properties of DC motors, 
where torque and motor magnetic flux can be independently controlled. In this 
scheme, the induction motor is analysed from a synchronously rotating reference 
frame where all the fundamental ac variables appear to be dc ones. The torque and 
flux component of currents are identified and controlled independently to achieve a 
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good dynamic response. However, there is a necessity of transforming the variables in 
the synchronously rotating reference frame to the stator reference frame to effect the 
control of actual currents/voltages. This transformation contains transcendental 
functions like sine, cosine, and also introduces computational complexity into the 
system. Additionally, the transformation also needs the appropriate flux vector angle, 
which is either calculated by adding the calculated slip angle and measured rotor 
angle as in indirect vector control [4, 5] or by estimating the flux angle directly by 
employing a flux observer as in direct vector control. Thus the accuracy of the vector 
control is largely governed by the accuracy with which flux angle is calculated and 
rotating reference frame variables are transformed into the stator variables. 

■ Field Acceleration method: This method is based on maintaining the amplitude of 
magnetic field in the air gap as constant and there is an adjustment in the speed of 
rotation of the field to produce the desired value of the torque. This also means when 
the air gap flux is kept constant, all other variables like voltage, current and torque 
become functions of slip frequency and torque control is achieved by accelerating and 
decelerating the rotation of flux [6, 7]. 

■ Direct Torque and Flux Control (DTFC) has emerged over the last decade to become 
one possible alternative to the well-known vector control of induction machines. Its 
main characteristic is the good performance, and it gives results as good as the 
classical vector control but with several advantages based on its simple structure and 
having four quadrant control capability. In DTFC [8-19] it is possible to control 
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directly the stator flux and the torque by selecting the appropriate inverter state. Thus 
the flux and torque are directly controlled by closed loop feedback as opposed to the 
vector control where the flux and torque are indirectly controlled by the direct and 
quadrature axis components of stator current vector. DTFC uses feedback control of 
torque and stator flux, which are computed with measured voltages and currents. The 
method uses stator reference frame model of the induction motor for its 
implementation, thereby avoiding trigonometric operations in the coordinate 
transformations of the synchronous reference frames [7]. The scheme uses stator flux 
linkages controls that enable the flux-weakening operation of the induction motor to 
be straight forward compared to rotor flux-weakening as the stator flux is directly 
proportional to the induced emf This reduces the dependence of the scheme on many 
motor parameters, thus making it robust, especially in the flux-weakening region. 
Thus the scheme depends only on stator resistance of the motor [20]. For flux and 
torque control, this scheme requires the position of the flux phasor, which is difficult 
to obtain at low and zero speeds from measured voltages and currents. At such low 
speed, the stator induced emfs are very small, making scaling and accurate 
problematic; also the stator-resistance variations introduce significant error in the 
flux-phasor computation. Thus the low-speed operation suffers with this scheme 
without stator resistance compensation [10, 15, 18]. There is also a possible 
instability problem, due to parameter sensitivity. The drive system becomes unstable 
if the measured stator resistance in the controller is higher than the actual machine 
resistance [22]. Although majority of applications are in the medium power range, the 
DTFC algorithm is equally applicable for the high-power motor drives. Most of the 
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three-level inverters are generally used in driving the high power induction motors. In 
the high power three level inverter applications, the inverter switching frequency is 
limited below 1 kHz [21]. This limitation on the switching frequency results from a 
consideration that the higher the switching frequency, the heavier the cooling 
apparatus for the respective switching semiconductors. As a result, such a 
consideration on the low switching frequency of the inverter is taken into account 
design the DTFC algorithm suitable for three-level inverter-fed induction motors [23, 
24, 25]. 

In the inverter-fed induction motor drives, a dc power source is required to feed 
the dc link. Traditionally, the power conversion from ac to dc was achieved by 
uncontrolled diode rectifiers or thyristorised controlled converters. Such converters have 
drawbacks such as harmonics in the input current, and output dc voltage, in addition to 
low input power factor for thyristors controlled converters particularly at low output 
voltages. Furthermore, four quadrant operation is not possible without employing a dual 
converter. Modern electric traction drives require the front-end ac to dc converter to have 
both rectifying and regenerating abilities. Dual converter used under such requirements 
has complexity in power and control circuits. The use of synchronous link three level 
converter (SLC) [26,27] as a front-end converter ensures near sinusoidal input current, 
reduced harmonic current injection into the utility and near unity power factor operation 
under forward and reverse power flow conditions. 
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1.2 Objectives of the Thesis 

The objective of the present work is to realize a variable speed three-level 
inverter- fed induction motor drive where the torque and the flux of an induction motor is 
controlled directly and independently by using Direct Torque and Flux control method. 
At the same time, unity input power factor is maintained by using a three phase three 
level SLC as the front-end ac to dc converter. The dc link voltage is regulated by a closed 
loop feedback using a PI controller. Unity power factor at the utility interface is obtained 
by controlling the front-end converter switches in such a way that the source currents 
follow the reference currents which are in phase with the respective source voltages. Thus 
the present drive is utility friendly. 

A PC is used to achieve the closed-loop operation. Thus, the present thesis deals 
with the design, simulation and PC-based implementation of a four quadrant direct torque 
and flux controlled induction motor drive system using three level inverter. 

1.3 Summary of the chapters 

The work presented in this thesis is organized in five chapters as given below: 

Chapter 2 describes the basic operation of a three phase three level synchronous 
link converter. Modeling and simulation has been made with resistive load. Simulation 
results are presented at the end of the chapter. 

Chapter 3 highlights the main feature of the proposed DTFC scheme along with 
description of its principle and basic concept of the torque and flux control. Mathematical 
modeling of induction motor and control algorithm of the proposed DTFC scheme using 
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three-level VSI is also presented in the chapter. Simulation has been made with the 
synchronous link converter used as the front-end converter. Typical simulation results are 
presented at the end of the chapter. 

Chapter 4 describes the hardware fabrication and PC-based implementation of 
control scheme using a high performance data acquisition card. Brief discussions about 
the power circuit and the control circuit are given with requisite diagrams. Experimental 
results are presented at the end of the chapter. 

Chapter 5 enumerates the main contributions of the present thesis and suggests the 
scope for the future work. 
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Chapter 2 


Modeling and Simulation of Three Level Synchronous Link Converter 


2.1 Introduction 

The power conversion from ac to dc has been traditionally dominated by 
uncontrolled rectifiers or phase controlled converters. Such converters have the inherent 
drawbacks such as, harmonics in the input current and output voltages, in addition to low 
input power factor particularly at low output voltages. Variable speed drive applications 
require ac to dc converters to have both rectifying and regenerating abilities. Dual 
converters can be used under such requirements, but they have complex power and 
control circuits. The Synchronous Link Converter (SLC) can be used as the front-end 
converter where sinusoidal input current with four quadrant operation is required. 

Synchronous Link Converter has several advantages such as sinusoidal input 
currents, unity power factor, inherent regenerative capability and these make SLC a 
plausible choice for ac to dc converter in a variable speed drive applications. In recent 
years, multilevel converters have become popular in medium and high power applications 
due to their ability to meet the increasing demand of power ratings. This chapter deals 
with the modeling and simulation of a three phase three-level Synchronous Link 
Converter, which can be used as a front-end converter for a variable speed drive system. 



2.2 Principle of Three-phase Three-level Synchronous Link Converter 

The power circuit configuration for the three-phase three-level is shown in Fig 
2.1. It is operated in controlled current boost type converter mode. The source currents 
) are made to follow sinusoidal reference currents within a hysteresis band. The 

hysteresis band width determines the source current profile, its harmonic spectrum and 
the switching frequency of the devices. From Fig. 2.1, the source currents 

be increased by turning on the upper two switches (A1 & A2, B1 & B2, Cl &C2) and 
source currents can be decreased by turning on the lower two switches (A3 & A4, B3 & 
B4, C3 &C4). The dc link capacitor voltage is kept constant throughout the operation. 



‘dc2 

Fig. 2. 1 Three-phase Three-level Synchronous Link Converter 
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The switching logic employing hysteresis current control can be written as follows. 

If zT > / r + A/ then switches A1 and A2, turn on. 

« a_rej 

If in <i r -AI then switches A3 and A4, turn on. 

« ajref 

Where AI is the hysteresis band on the reference current. Similarly, the same logic 
can be applied for other phase currents also. 

2.3 Hysteresis Current Control Scheme 



Fig 2.2 Hysteresis current control scheme for SLC 
The hysteresis current control scheme block diagram is shown in Fig. 2.2. The dc 
link voltage is compared with the reference dc link voltage to control the amplitude of the 
reference ■ currents. The error between the actual dc link voltage and reference dc link 
voltage is processed through a PI controller, which gives the amplitude of the reference 
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currents. The amplitude is multiplied by sin(wt), sin(wt-120), sin(wt-240) respectively, 
which are in phase with the respective phase voltages, to generate the reference source 
currents of the converter. The switches are controlled in such a way that the actual 
currents follow the reference currents with in the present hysteresis band. 

2.4 Modeling of Three-phase Three-level SLC 

The modeling of SLC is divided into two stages. In the first stage, we obtain the 
expressions for the current dynamics, and in the second, we obtain the expressions to 
describe the dynamics of the capacitor voltages. 

A. The inductor current dynamics 
Direct application of Kirchoff s voltage law (KVL) yields, 
di. 


'^as - ^s + ^r^a + '*^<30 + ^on 


(2.1) 


"6s 


dih 


( 2 . 2 ) 


dir 

^cs = + ^r^c + ^co + ^on 


(2.3) 


Since source voltages are balanced and a three phase three wire system is 

considered, we can write the following equations 

V + V, + V =0 
as bs cs 


i + iT + i = 0 
a b c 


Adding all the three equations (2.1), (2.2) and (2.3), we get 


^on = --(v +v, +V ) 

3 bo CO 


Eliminating from (2.1), (2.2) and (2.3), we get 
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+ R ii 
r b 


■1 -1 


■1 -1 2 


The pole voltages of the SLC (where i = a, b, c) with respect to point 0 are 
defined in terms of the switching variables S.-^ and S.^ as follows: 


V. = 5 V ^ ^ 

10 \ fi dc\ i2 del 


[1 when switches A: and^ are on 
5 = 1 2 

fi [ 0 else 

fl when switches A: dsidk , are on 
S.^=\ 3 4 

0 else 


where k = A,B,C and i = a,b,c 


Fi'om equation (2.4) we can get, 


•1 -iTv, 


^bs ib “T ^ 


' ^ ''bo 

1-1 2 V(<Q 


Hence, the inductor current dynamics can be modeled as (2.6) where the SLC 
pole voltages v^q , and are computed in (2.5). 

B. The capacitor voltage dynamics: 

The dc link currents i^^^ and i ^^2 evaluated using Kirchoff s Current Law 


(KCL) as follows. 
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^dc\ ^ ^ k\^ k 

k = a 


^dc2 , ^ ^ kl\ 

k = a 


’k\ 




1 when switches A:, and k 2 are on 
0 else 


1 when switches k and k. are on 


0 


else 


C, 


dv 


dc\ 


dt 


= i 


dc\ 


^dc\ 

R, 


and C. 


dv 


del _ ■ _ ^dc2 

~ '■del 


dt 


R, 


(2.7) 


( 2 . 8 ) 


(2.9) 


Ci = C2=C 

We can simplify equation (2.9) as, 


dv 


dc\ 

1 

dt 

dVdc2 

C 

_ dt _ 



'rfcl 


^dc2 


''dc\ 

Ri 

^dc2 

Ri 


( 2 . 10 ) 


Hence, the capacitor voltage dynamics can be modeled as (2.10) where the current 
and given in equations (2.7) and (2.8) respectively. 

The equations (2.6) and (2.10) can be solved by R-K fourth order method and the 
values of , Vdc2JaJb and can be obtained. is the sum of and Vdc2- 

A computer program has been developed for resistive load. The data for 
simulation are given below. 

1 . Input supply voltage = 230 V (rms), 30, 50 Hz. 

2. Output dc voltage = 350 V. 

3 . Synchronous Link Inductor =8 mH. 

4. DC link Capacitors C, =0^ =1000/zF. 

5. Load resistance i?/ = 20 Q. 
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In the simulation results, the input phase current waveforms, the source voltages, 
and the dc link voltages are shown in following figures. Fig. 2.3 shows that the dc link 
voltage is nearly constant during the operation. Fig 2.4 and 2.5 show that the sharing of 
voltage between capacitors is same. Figs. 2.6, 2.7 and 2.8 show that input phase currents 
are in phase with the respective voltages. 

2.5 Conclusion 

Due to advantages like reduced line current harmonics, unity power factor 
operation and bi-directional power flow capability, synchronous link converters are 
widely used as front-end converters in regenerative ac drives. Hysteresis current control 
scheme is used, and unity power factor operation is shown in simulation results. Usage 
of multilevel converter has increased the power rating and reduced the harmonic 
distortion. The use of multilevel converter in the front-end allows transformer less 
interface with high voltage utility for high power application. 

2.6 Simulation results : 



Y-axis Idiv = 100 V ( DC voltage) 
Fig. 2.3 Response of DC link voltage 
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Fig. 2.4 Response of DC link voltage 



Fig. 2.5 Response of DC link voltage 
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Vas,la vs Time 



X-axis Idiv = 0.05 sec (time) 

Y-axis Idiv = 100 V ( source voltage, ) 

Y-axis Idiv = 10 A (source current, i^) 

Fig. 2.6 Waveform of source voltage and source current for phase ‘a’ 



y-axis Idiv = 100 V ( source voltage, ) 


y-axis Idiv = 10 A (source current, ) 

Fig. 2.7 Waveform of source voltage and source current for phase ‘b’ 
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iajb.k: (A) 



X-axis Idiv = 0.05 sec (time) 
Y-axis Idiv = 100 V ( source voltage, 


Y-axis Idiv = 10 A (source current, i^) 

Fig. 2.8 Waveform of source voltage and source current for phase ‘c’ 



X-axis Idiv = 0.05 sec (time) 

Y -axis 1 div = 1 0 A (source currents, ) 

Fig. 2.9 Waveform of source currents 
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Chapter 3 


Modeling and Simulation of Converter-Inverted fed Induction Motor 
Drive with Direct Torque and Flux Control 


3.1 Introduction 

The various control strategies for the control of the inverter-fed induction motor 
have resulted in good steady-state but poor dynamic response. From the plots of the 
dynamic responses, the cause of such poor dynamic response is found to be the deviation 
of the air gap flux linkages from their set values. The deviation is not only in magnitude 
but also in phase. The variations in the flux linkages have to be controlled by the 
magnitude and frequency of the stator and rotor phase currents and their instantaneous 
phases. Scalar control strategies have utilized the stator phase current magnitude and 
frequency and not their phases. This resulted in the deviation of the phase and 
magnitudes of the air gap flux linkages from their set values. 

The oscillations in the air gap flux linkages results in oscillations in 
electromagnetic torque and, if left unchecked, result as speed oscillations. This is 
undesirable in many high-performance applications, such as in robotic actuators, 
centrifuges, servos, process drives, and rolling mills, where high precision, fast 
positioning, or speed control is required. Such requirement will not be met with the 
sluggishness of control due to the flux oscillations. Further, air gap flux variations 




resulted in large excursions of stator currents, requiring large peak converter and inverter 
ratings to meet the dynamics. An enhancement of peak inverter rating increases cost and 
reduces the competitive edge of ac drives in marketplace, in spite of excellent advantages 
of the ac drives over dc drives. 

Separately excited dc drives are simpler in control because they independently 
control flux, which, when maintained constant, contributes to an independent control of 
torque. This is made possible with separate control of field and armature current 
magnitudes, providing simplicity not possible with ac machine control. In contrast, ac 
induction motor drives require a coordinated control of stator current magnitudes, 
frequencies and their phases, making it a complex control. As with the dc drives, 
independent control of the flux and torque is possible in ac drives. The stator current 
phasor can be resolved, say, along the rotor flux linkages, and the component along the 
rotor flux linkages in the field-producing component, but this requires the position of the 
rotor flux linkages at every instant as it is dynamic unlike the dc machine. The 
requirement of phase, frequency, and magnitude control of the currents and hence of the 
flux phasor is made possible by inverter control. 

Vector control made the ac drives equivalent to dc drives in the independent 
control of flux and torque and superior to them in their dynamic performance, because it 
relates to the phasor control of the rotor flux linkages. In vector control method, the 
machine control is achieved in a synchronously rotating reference frame in which the ac 
variables appear as dc quantities. The torque and flux components of stator current are 
identified and controlled independently to achieve good dynamic response. However in 
vector control there is a necessity of transforming the variables in synchronously rotating 
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reference frame to stator reference frame, to effect the control of actual motor 
current/voltages. This transformation requires trigonometric functions and so introduces 
computational complexity into the control system. Additionally, the transformation also 
needs the flux vector angle which is either calculated by adding the calculated slip angle 
and measured rotor angle as in the indirect control or by estimating the flux directed by 
employing a flux observer. Thus, the accuracy of vector control id largely governed by 
the accuracy with which the flux angle is calculated and rotating frame variables are 
transformed into the stator variables. 

In recent years, the commercial applications of the field- oriented induction motor 
drives have greatly increased. DTFC strategy for ac motor control seems to be 
particularly useful for the adjustable speed drives, due to its robustness and functional 
simplicity. Despite its simplicity, DTFC is able to produce very fast torque and flux 
control and is robust with respect to motor parameters changes and to perturbations. 
Unlike the conventional vector control, DTFC requires no coordinate transformation of 
variables. Using the space-vector modulation (SVM) provides higher control resolution 
and helps improving the drive’s behavior. The DTFC for ac drives is a strategy 
exclusively based on stator voltage control. The consecutive voltage vectors applied to 
the motor are directly selected on the basis of torque and flux errors. In this way, fast 
response and robust torque and flux control are obtained, without intermediate current 
control. The classic DTFC uses bang-bang torque and flux controllers, without 
decoupling. A simple switching logic (switching table) employs the output signals of 
these controllers to select the most appropriate voltage vector, i.e., the one which rapidly 
reduces the torque and flux errors. 
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High power induction motor drives using classical three-phase converters have 
the disadvantages of poor voltage and current qualities. To improve these values, the 
switching frequency has to be raised which causes additional switching losses. Another 
possibility is to put a motor input filter between the converter and motor, which causes 
additional weight. A further inconvenience is the limited voltage that can be applied to 
the induction motor determined by the blocking voltage of the semiconductor switches. 
For high power semiconductors, the switching frequency is limited by the maximal 
power loss. The aim of three-level converter is to generate motor currents of high quality 
by using more semi conducting devices. In particular, the converter is used for high- 
voltage motor drives. The series connection allows reaching much higher voltages than 
the blocking voltage of the semiconductors. The resulting switching frequency is a 
multiple of the switching frequency applied to the switches. Variable speed induction 
motors need regenerative braking capability. If diode or phase controlled rectifiers are 
used as a front end converter, regenerative braking is not possible. In addition, for these 
converters source current has harmonics, and the power factor is less then unity. The 
above problem can be solved by using the synchronous link three-level converter as the 
front-end converter. 

The present chapter is devoted to the development of control scheme and digital 
simulation of direct torque and flux control of induction motor using three-level inverter 
with current control synchronous link converter as front-end converter. The induction 
motor is simulated in the stationary reference frame. Typical simulation results are given 
at the end of the chapter. 
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3.2 Converter - Inverter fed Induction Motor Drive 

Fig. 3.1 shows the basic block diagram of three-level converter-inverter fed 
induction motor drive system. The drive system consists of two major parts: the driving 
part that contains the ac to dc three-level converter, three-level dc to ac inverter, and the 
induction motor, as well as the control part consists of the controller circuit and software. 
The three-level ac to dc converter is called synchronous link converter (SLC) to make the 
system utility friendly as discussed in chapter 2. Three-level inverter is used to control 
the induction motor. In the following sessions details of direct torque control of 
induction motor using three-level inverter will be given. Initially the dynamic modeling 
of induction motor has been done. 


AC supply 
30 



Fig 3.1 Block diagram of three-level converter-inverter fed induction motor drive 
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3.3 Dynamic Modeling of Induction Motor 

A dynamic model of the machine subjected to control must be known in order to 
understand and design controlled drives. Due to the fact that every good control has to 
face any possible change of the plant, it could be said that the dynamic model of the 
machine could be just a good approximation of the real plant. Nevertheless, the model 
should incorporate all the important dynamic effects occurring during both steady-state 
and transient operations. Furthermore, it should be valid for any changes in the inverter’s 
supply such as voltages or currents. 

The dynamic model of the induction motor is derived by using a two-phase motor 
in direct and quadrature axes. The concept of power invariance is utilized in the 
modeling. According to the concept of power invariance, the power must be equal in the 
three-phase machine and its two-phase equivalent model. 

The following assumptions are made to derive the dynamic model: 

• Uniform air gap; 

• Balanced rotor and stator windings with sinusoidally distributed mmf; 

• Inductance vs. rotor position is sinusoidal; and 

• Saturation and parameter changes are neglected. 
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3.3.1 Modeling of Induction Motor in Arbitrary Reference Frame 

The dynamic equations of the induction motor in arbitrary reference frame can be 
represented by using flux linkages as variables. The stator and rotor flux linkages m 

arbitrary reference frame are defined by 

fC _ T jC r :C (3.1) 

fds " ^ds ^ dr 

rC r I r (3.2) 

fqs = ^s^qs + ^m^qr 

fC _ r jC , r .-C (3.3) 

fdr - 'dr + 'ds 

rC ^ r ^ T f (3-4' 

Jqr ~ ^r ^qr "*■ ^qs 

The stator and rotor voltages can be written as 


11 

Vi ^ 

^ds ~ 

+ Pfds 

(3.5) 

-- 

^ qs 

= ^qs Rs f ds 

+ P fqs 

(3.6) 

II 


\fqr + Pfdr 

(3.7) 

vqr - 


)^dr 

(3.8) 


By substituting equations 3.1 - 3.4 into equations 3.5 - 3.8, the induction motor 
model in arbitrary reference frame is obtained. For simplicity superscript “c” is removed 
in the matrix. 

"^qs ^s'^^sP ^n^P 

-OcLs Rs+^sP ^mP ^ds (39) 

0^ " LfnP [cOc-cOrYm Rr + P^r (^c-^r)Rr V 
0 "r )Rm ^mP ~ i^c “ ^r )^r Rr + P^r Ijdr . 
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3.3.2 Modeling of Induction Motor in Stationary Reference Frame 

The motor is fed from voltage source inverter as shown in fig. (3.2). For a given 

set of switching information, V , V, and V have discrete instantaneous values. 

° CIO ho CO 


Thus, 


Van = --(2*1^ -F, -V ) 

3 ^ ao bo CO 


(3.10) 


V, = --(2*F, -V -V ) 

[jyi ■a'* 


ho ao CO ‘ 


(3.11) 


Vcn = --{2*V -V, -V ) 


(3.12) 


3 ' CO ho ao 
Assuming that the direct phase winding axis (d-axis) of two-phase model is along 
the a-phase winding axis of the three-phase machine, transformation of ac voltage and 
current from a, b, c (3 -phase) frame to d, q (2-phase) stationary reference frame gives 

■ ^ds-^Ahn-oAhn^^cn)) (3-13) 

\V^..-Vrn) 


Vqs 


Aj3y ^ 


(3.14) 

(3.15) 


^qs 




(3.16) 


The DTFC scheme uses a stator reference frame model of the induction motor for 
its implementation. Therefore stator reference frame model is obtained by putting co^ = 0 
in the arbitrary reference model of induction motor. 
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+ 

=0 

11 

(3.18) 


^qs ~ ^qs P fqs 

(3.19) 

We can get 


(3.20) 


fgs ~ 1 ^gs ~ ^gs^s )‘2t 

(3.21) 


Therefore, 


( 3 - 22 ) 

The position of stator flux and hence the sector can be obtained as 



The electromagnetic torque can be expressed as: 

=(V2lPI2lf^igs-fisia,) ( 3 . 24 ) 

The speed of the rotor in electrical rad/sec is 

( 0 , = (P/ 2 X 1 /JX 3 ; - 2) - BiZ/PK ) (3.25) 
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3.4 Direct Torque and Flux Control 


DTFC main features are as follows; 

• Direct control of flux and torque, which are computed from measured 
values of stator voltages and currents. 

• Indirect control of stator currents and voltages. 

• Approximately sinusoidal stator fluxes and stator currents. 

• High dynamic performance. 

The main advantages of DTFC are: 

• Absence of coordinate transformation. 

• Absences of voltage modulator block, as well as other controllers like PID 
for motor flux and torque. 

• Minimal torque response time even better than vector controllers. 

However, some disadvantages are also present such as: 

• Possible problems during starting. 

• Requirement of flux and torque estimators, implying the consequent 
parameters identification. 

• Inherent torque and stator flux ripple. 
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3.4.1 Basic Concepts of DTFC 


The torque expression given in Equation (3.24) can be expressed in the vector 


form as 


T =- 
2 


fp\ 






(3.27) 


where fs=f ~ J ~ J this equation, is to be replaced by 

S CIS S CIS 

rotor flux fy . In the complex form, and fy can be expressed as functions of currents 


as 


fo=L 1 +L I 
^ s s m r 


(3.28) 


fy^L I +L 1 


r r ms 


(3.29) 


Eliminating ly from Equation (3.28), we get 


- L - , - 

f =-!lLf +L I 

JS Jy s s 

r 


(3.30) 


where ~ rn * corresponding expression of is 


, _ 1 r r 

- tJs TJr 


(3.31) 


Substituting Equations (3.31) in (3.27) and simplifying yields 


^ 2l 2 


■^m 


LyLg 


'fr'^fs 


(3.32) 


that is, the magnitude of the torque is 
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phasor (or vector) diagram for Equation (3.32), indicating the vectors , fy and Ig for 
positive developed torque. If the rotor flux remains constant and the stator flux is 


changed incrementally by stator voltage Fy as shown, and the corresponding change of y 
angle is A/ , the incremental torque NTq expression is given as 


P^ L, 


Ly Ay 


\fr\ fs + ^fs\ 


(3.34) 



d-axis 

Fig. 3.2 stator flux, rotor flux, and stator current vectors on d-q stationary 
reference plane (stator resistance neglected) 
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3.5 Three-Level Space Vector Modulation 

Multilevel voltage-fed inverters with space vector modulation have established 
their importance in high power and high performance industrial drive applications. Three- 
level voltage fed inverters are recently showing growing popularity for multi-megawatt 
industrial drive applications. The main reasons for this popularity are easy sharing of 
large voltage between the series devices and improvement of harmonic quality at the 
output compared to a two-level inverter. The space vector modulation of a three-level 
inverter provides the additional advantage of superior harmonic quality. 



Fig. 3.3 Three-Level Inverter fed Induction Motor 
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However, space vector modulation of a three-level inverter is considerably more 
complex than that of a two-level inverter because of large number of inverter switching 
states. In addition, there is a problem of neutral point voltage balancing. 

The structure of a three-level neutral point clamped inverter is shown in Fig. 3.3. 
A three-level inverter differs from a conventional two-level inverter in the sense that it is 
capable of producing three different levels of output phase voltage. When switches 1 and 
2 are on, the output is connected to the positive supply bus. When switches 3 and 4 are 
on, the output is connected to the negative supply bus. When switches 2 and 3 are on, the 
output is connected to the supply neutral point via one of the two clamping diodes. With 
three possible output states for each of the three phases, there are a total of 27 possible 
switch combinations. 



Fig. 3.4 Switching states of three-level inverter 
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The result of plotting each of the output voltages in a d-q reference frame is 
shown ill Figure 3.4. Figure 3.4 shows that the 27 switch combinations result in a total of 
19 unique voltage vectors since some of the combinations produce the same voltage 
vector. These different combinations relate to different ways of connecting the load to the 
DC bus that result in the same voltage being applied to the motor. 

3.6 Theory of DTFC Involving Three - Level Inverter 

In this section, the basic principle for the direct torque and flux control in a three- 
level inverter system is explained. Fig. 3.5 shows space-vector representation of three- 
level inverter output voltages, in which subscripts z, h, i and f denote zero, half, 
intermediate, and full voltages, respectively. 

I 



Fig. 3.5 Voltage vectors of three-level inverter 
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Since 27 possible choices for switching voltage vector selection exist in the three- 
level inverters, appropriate selection of the inverter switching voltage vector is more 
complex than that in the two-level inverters. Assuming that the stator flux' vector is in the 
1^‘ sector, selection of the respective stator voltage vector is described in Fig. 3.5. The 
selection of ^2 f h ^2 f h increase the phase angle between the 

corresponding stator and rotor fluxes. As a consequence, the developed torque can be 
increased by the application of these voltages. It can be seen in Fig. 3.5 that the stator 
flux is increased by the selection of, ^2 f h decreased by y ^ . If the half voltage 

vector is selected, the lower slope of torque can be obtained. 

Upper 
torque band 
Lower 
torque band 

Torque Reference 

- Lower 
torque band 

- Upper 
torque band 

► 2- Level torque slope 

► 3- Level torque slope 

Fig. 3.6 Torque slope pattern of three-level inverter 

Fig. 3.6 shows a double torque hysteresis band method applied for three-level 

inverter. It is seen in this figure that the torque control method for three-level inverter is 

provided with the lower torque bands additionally to the hysteresis bands for two-level 

inverters. 
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fs in sector k 

Torque 

t 

tt 

i 

Flux 

t 

^k+l,h 

'^k+fif 


i 

'^k+2,h 

^k+2,f 

^z 


Table 3.1 

Fig. 3.7 Switching lookup table 


It is noted that positive lower hysteresis band is used for forward rotation of the 
induction motors and negative lower band is for reverse rotation. Provided that the 
induction motor is rotating in the forward direction, for instance, the appropriate full 
voltage vector is chosen to increase the developed torque when instantaneous torque goes 
down to the negative upper hysteresis band. When the controlled torque reaches the 
positive lower hysteresis band, the full voltage vector is replaced with the respective half 
voltage. If torque increases beyond the positive upper torque band, the zero voltage 
vector is applied to decrease the developed torque value. For the reverse rotation, 
considering that the reverse rotating field is required, the rules similar to the selection for 
the forward rotation can be applied. The resulting switching lookup table is shown in 
Figure 3.7. 

When the above-mentioned DTFC switching method is applied for the low speed 
operation, it causes the demagnetization phenomenon of the stator flux, as shown in Fig. 
3.8. This problem can be explained as follows. The demagnetization by nonzero-voltage 
vector is described in Fig. 3.8 (a), where the stator flux is placed in the k+l^"^ sector. It is 
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seen in Fig. 3.8 (a) that, f^, the stator flux in the region near the sector-to-sector 
boundary, happens to be decreased by the stator voltage vector ^”^^2 f h 

selected to increase the stator flux magnitude by the basic switching lookup table of Fig. 
3.7. It is noted that the voltage drop across the stator resistance takes significant 
proportion of the applied voltage, especially in the low-speed region. 

3.7 Switching Strategy for Low - Speed Operation 

The basic DTFC algorithm has some problems such as the stator flux drooping 
and the poor torque control, especially in the low-speed region. The flux drooping 
problem appears just after the position of the stator flux vector moves from one sector 
into the next sector. The application of the zero-voltage vector makes it difficult to 
decrease the developed torque in the low-speed operation. 

The basic switching look up table described in the previous section cannot solve 
these problems occurring in the low-speed operation. Even if the flux demagnetization 
problems can be worked out by a rotation of the reference frame or the switching sector 
that are similar to methods for two-level inverters, a new method using the intermediate 
voltage vectors is more desirable for the three-level inverters. Moreover, the intermediate 
voltage vectors can be effective for preventing the torque control from being inactive. 
Fig. 3.9 shows only one sector of the 12 sectors that result from subdividing a six-sector 
plane of Fig. 3.5 into 12-sector plane. The every sector in Fig. 3.5 subdivided into the 
lower sub sector and the upper sub sector, each of which has a width of 30°. It is assumed 
that the stator flux rotates in the counter clockwise direction. 
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(a) Demagnetization by nonzero voltage vectors 



(b) Demagnetization by zero voltage vectors 
Fig. 3.8 Demagnetization problem in low-speed operation 
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Fig. 3.9 Improved switching method 


In the lower sub sector, the stator voltage vector y ^ which is determined 

according to table Fig. 3.7, becomes ineffective in boosting the stator flux because the 
full voltage vector is not able to sufficiently provide a flux vector component to increase 
the corresponding stator flux. It can be seen in Fig. 3.9 that the respective intermediate 
stator voltage . can be more effective to increase the resultant stator flux. If the 


intermediate voltage vector . is selected instead of F^_^j ^ , the demagnetization 

problem can be resolved. When the stator flux moves into the upper sub sector passing 
the lower sub sector, the voltage selected, F^^^ f h’ increase the developed 

torque and the voltage vector F^_^ y ^ instead of the zero voltage vector being used to 

decrease the developed torque. 
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Torque band 



(a) Switching pattern for low half sector 


Upper 
torque band 
Lower 
torque band 

Torque reference 


-Lower 
torque band 
-Upper 
torque band 



(b) Switching pattern for upper half sector 
Fig. 3.10 Torque slope pattern in low-speed operation 


fg in sector k 

Torque 

T 

n 

i 

it 

Flux 

Upper 

subsector 

^k+l,h 

T)k+i,f 

T>k-l,h 

Dk-i,f 

Lower 

subsector 

^k,i 

Dk-2,i 


Table 3.2 

Fig. 3.11 Modified lookup table 


To obtain the effectively magnetizing effect around the boundary between two 
sectors, especially at the low-speed region, a modified lookup table as shown in Fig. 3.11 
is devised. The torque slope pattern change is shown in Fig. 3.10, where two and four 
torque slopes are applied in the lower subsector and the upper subsector, respectively. 
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3.8 Application of Synchronous Link Converter in DTFC of Induction Motor 

The analysis and simulation of the synchronous link converter with resistive load 
has been discussed in chapter 2. When the synchronous link converter used as the front 
end converter in induction motor drive system, the power balance equation of the system 
is given as follows. 

inverter input = motor input ( by neglecting converter losses) 


rotor input + stator losses 


V, */ =1 *0) +(3/2)*/ *I *R 

dc ac e /- w / ^ 


(3.35) 


where <x>^ is the motor speed, and is electromagnetic torque developed by 


induction motor as given by Equation (3.24). 

The load current for the synchronous link converter is given by 

I , =[t *co +(3/2)*/ */ *R )/f, 

dc ^ e r ^ ‘ ’ rms rms sij dc 


(3.36) 


In the motoring mode the inverter input current f^^is positive, i.e. the power 
flow is from the source to the motor and the input power factor is maintained at unity as 


shown in simulation results. In the regenerative mode current / , is negative i.e. the 


induction motor act as a generator and the power is fed from the motor to the source. 
During this period source current and source voltage are out of phase (180°). 
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3.9 Four Quadrant DTFC Operation of Induction motor 

To visualize the different steps of the simulation, a block diagram has been drawn 
as shown in Fig. 3.13. The source current reference is obtained by the dc link voltage 
error processed through a PI controller as given as 




The reference torque is generated by the speed error processed through 


the PI controller. The output of the PI controller is given as follows 


T .W = r An l)+K (e -e / + r ,\A/ (3.38) 

e_rej^ ' ~ ’ p\ wrn wr\n-\)) i wr\n-\) 

The model block takes the three-phase motor currents, the switching status of the 
inverter, the dc link voltage, and calculates the stator flux and the torque. The actual 

values of the electromagnetic torque Tq and actual flux / are compared with the 

s 

reference torque T ^ r and reference flux/ r to give the torque and flux errors. 

(i T Sj S 

Now the torque and flux errors are digitized by the five-level and two-level hysteresis 
comparators respectively. The outputs of the 5-level hysteresis torque comparator ct, 2- 

level flux comparator cf, and sector angle command 9^ is given to the optimum switching 

table which generates the required switching pulses for the inverter as the table shown in 
the Fig. 3.12. 
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Fig. 3.13 Block diagram of DTFC with SLC as front end converter 
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Induction Motor 












3.10 Simulation Results 

The flow chart of the computer simulation is shown in Fig.3.14, and the program 
listing in C language is given in the Appendix D. The inputs to the program are motor 
parameters like stator and rotor resistances, inductances, moment of inertia of the coupled 
drive system, reference flux, reference speed, and reference dc link voltage. The 
simulation plots are given in next few pages. 

Fig 3.15 (a) shows the demagnetization of the stator flux under low speed 
operation (5% of base speed) by using conventional look up table. It can also be 
observed in Fig. 3.16 (a), polar plot of dq-axes stator flux. As from the Fig. 3.15 (b) and 
3.16 (b) we can conclude that the flux drooping is almost removed and stator flux follows 
a circular trajectory even in the low speed region. 

Fig. 3.17 to Fig. 3.27 show the various wave forms under no load operation at 
85% of the rated speed of the motor. Fig. 3.20 shows the dc link voltage obtained by 
using three-level synchronous link converter. Fig. 3.21 shows that unity power has been 
achieved by using synchronous link converter and Fig. 3.22 shows that the source 
currents are sinusoidal. 

Fig. 3.28 to Fig. 3.39 shows the motor currents, flux, speed and torque responses 
at the starting on no load, and various responses to step change in speed from standstill to 
750 rpm at starting, speed reversal from 750 rpm to - 750 rpm at 1.75 second and again 
the speed reversal from reverse motoring (-750 rpm) to forward motoring (750 rpm) at 3 
second. Speed reversal from forward motoring at 750 rpm to reverse motoring -750 rpm 
takes around 0.25 second. Motor current is found to be in limits during starting and 
during speed reversal. Motor currents and d, q-axes fluxes are nearly sinusoidal. From 
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Fig. 3.36 it can be seen that there is a drop in the DC link voltage during transition of 
speed. It is clear from Fig. 3.37 that unity power factor has been achieved under forward 
motoring operation. Fig 3.38 shows that unity power factor is still maintained under 
reverse motoring operation. It can be observed that the power generated under reverse 
motoring operation is dissipated in the SLC. 

Fig. 3.40 to Fig. 3.50 shows motor currents, flux, speed and torque responses on 
application and removal of a load of 2 N-m at 1 .5 second and 4 second respectively. 
Torque response is found to be almost instantaneous. Speed of the motor follows the 
reference speed with a dip and a raise corresponding to change in load. Finally speed 
settles at the reference speed. Motor current is found to be in limits on application and 
removal of a load torque also motor currents and dq-axes fluxes are nearly sinusoidal. 
Fig. 3.48 shows that there is drop and raise in DC link voltage at the instants of change in 
torque. It can be observed that unity power is maintained under load condition from Fig. 
3.49. The source currents are being maintained approximately sinusoidal. No change in 
flux is observed with change in load and therefore independent control of torque and flux 
is achieved. 

3.11 Conclusion 

The analysis of direct torque and flux control of an induction motor has been 
made in this chapter. Synchronous link converter is used as the front-end converter for bi- 
directional power flow and to get unity power factor at the source side. The machine 
parameters are assumed to be constant at all operating conditions. The torque response 
depends upon the time differential of the slip angle and upon the capability of voltage 
source inverter to meet the required large transient power. Hysteresis current controlled 
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scheme is used for synchronous link converter. An advanced DTFC algorithm with 
improved low-speed operation of induction motors for high-power three-level 
applications has been presented in this chapter. The simulation results of induction motor 
drive under direct torque and flux control show excellent results in both steady state and 
transient conditions. 



Fig. 3.14 Flow chart for the Four Quadrant operation of DTFC 
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Flux.fs (Wb) 



(a) Using conventional lookup table 



(b) Using modified lookup table 


Fig. 3.15 Stator flux under low speed (5% of rated speed) operation 
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(a) Using Conventional lookup table 



(b) Using modified look up table 

Fig. 3.16 Polar plot of d-axis and q-axis stator flux for 5% of rated speed operation 
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Actual speed (Series 1), Reference speed (Series 2) 
Fig.3. 17 No load Operation of the motor with 85% of rated speed 



Fig. 3.18 Torque under no load at 85% of rated speed 
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Fig. 3.19 Stator Flux under no load and at 85 % of rated speed 
(Normal operating condition) 



Reference voltage (Series 1), Actual voltage (Series 2) 
Fig. 3.20 DC link Voltage at 85 % of rated speed operation 
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X-axis Idiv = 0.02 sec (time) 

Y-axis Idiv = 20 A (source current, ) — (Series 1) 


Y-axis Idiv = 100 V ( source voltage, v ) — (Series 2) 

CIS 


Fig. 3.21 Source Voltage and Source current at 85% of rated speed operation of the motor 



Fig 3.22 Source phase currents i „ under steady state operation 

CZtj Uu OS 
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Fig. 3.25 q-axis Stator flux under normal operating condition 



Fig. 3.26 d-axis Stator flux under normal operating condition 
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Speed (rpm 



Reference speed (Series 1), Actual speed (Series 2) 

Fig. 3.28 Forward motoring and reverse motoring Operation of DTFC 



Fig. 3.29 Stator flux during speed reversal 
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Fig. 3.30 Torque Response under speed changes 



Fig. 3. 3 1 d-axis or a-phase motor current during speed changes 
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fds (Wb) 



Fig,3.33 d-axis stator flux during speed changes 
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Reference voltage (Series 1), Actual voltage (Series 2) 
Fig. 3.36 DC link Voltage under step change in speed 



Y-axis Idiv = 20 A (source current ) 


Y-axis Idiv = 100 V ( source voltage, 

Fig. 3.37 Source Voltage and Source current (phase ‘b’) at 750 rpm (forward motoring ) 
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Y-axis Idiv - 20 A (source current i ) 

cs 


Y-axis Idiv = 100 V ( source voltage, v ) 

OS 

Fig. 3.38 Source Voltage and Source current (phase ‘c’) at -750 rpm (reverse motoring) 
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Actual torque (Series 1), Reference torque (Series 2) 
Fig.3.4() Step change in load torque at 705 rpm speed 



Fig.3.41 Speed during step change in load torque 
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Fig. 3. 42 Stator flux during step change in load torque 



Fig. 3.43 d-axis stator current during change in load torque 
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Fig. 3.44 q-axis stator current during change in load torque 



Fig. 3.45 d-axis stator flux during change in load torque 
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Fig. 3.46 q-axis stator flux during change in load torque 



Fig.3.47 Polar plot of d and q axis stator flux under change in load torque 
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Actual voltage (Series 2), Reference voltage (Series 1) 
Fig. 3. 48 DC link voltage during step change in load torque 



Y-axis Idiv = 20 A (source current, i ) 

. Cltj 

Y-axis Idiv = 100 V ( source voltage, 

Fig. 3.49 Source Voltage and Source current (phase ‘b’) under loaded condition 
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Fig 3.50 Source phase currents under loading condition 

(Steady state) 
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Chapter 4 


PC-based Implementation of the Induction Motor Drive System 


4.1 Introduction 

This chapter describes the experimental setup of PC-based hardware 
implementation of the “Direct Torque and Flux Control” of induction motor using three- 
level VSl, as discussed in chapter 3. In chapter 3, simulation studies have been made for 
DTFC of induction motor drive system. After obtaining the simulation results, it is 
important to validate the control algorithm under real time applications. This chapter 
gives a detailed description of experimental setup. Section 4.2 gives a detailed 
description of the experimental setup and the circuits used. Section 4.3, describes the PC 
- based implementation of the control circuit and section 4.4, PC - based control of the 
proposed DTFC scheme. Experimental results are given at the end of the chapter while 
section 4.6, ends this chapter with the necessary conclusions. 

4.2 Experimental Set - up 

The lay out of the PC - based implementation is given in Fig. 4.1. The practical 
implantation is carried out in such a way that the hardware circuits are minimized. 

The DC link voltage is supplied with a 3- phase rectifier. The 3-phase rectifier is 
fed from a 400 V, 3 phase AC supply through an auto transformer. Two 2200 pF 



capacitors are connected in series across output of the rectifier. Mid point of the two 
capacitors is treated as neutral. 

A motor generator set is used for the experiment. The detailed rating of the motor 
generator set is given in Appendix A. The parameters of the- machine are calculated by 
no-load and blocked rotor tests. Moment of inertia is found out from retardation test. 



j, 

Fig. 4. 1 Schematic diagram of the PC -{^sed implementation 

i 

For inverter switches, Mitsubishi IGBT modules (2 IGBTs in each module) with a 
voltage rating of 1200 V and current rating of 50 A, have been used. A high performance 
Mitsubishi driver circuit is used to drive IGBT. Actually it is a hybrid IC (M57959L) 
with a circuit designed for driving n-channel IGBT modules in any gate amplifier 
application. This device operates as an isolation amplifier for these modules and provides 
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the required electrical isolation between the input and output with an opto -coupler. Short 
circuit protection is provided by a built in de-saturation detector. The pin diagram and 
connection details for these drivers are given in the table 4.1 and Figures 4.2, and 4.3 
respectively. 

Table 4. 1 : Pin Configuration of the Hybrid IC M57959L 
Pin 1 : connected to collector of IGBT through a diode (1N4007) to protect against 
malfunction of over current protection (Detect pin) 

Pin 2: Not connected 
Pin 3: Not connected 
Pin 4: + 15 V supply (Vcc) 

Pin 5: Drive output given to the Gate of the IGBT 
Pin 6: - 9 V supply (Vi.:,.;) 

Pin 7-12: Not connected 
Pin 13: Drive Signal i/p (-) 

Pin 14: Drive Signal i/p (+) 


Block Otagram 



Fig. 4.2 Functional Block Diagram of IC M57959L 
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It Cdii bt nottd that to build a three-level inverter 12 IGBTs are required. So there 
is a need of 12 IGB'F driver circuits. Each Driver circuit requires +15 V and - 9 V dc 
voltage.s. 12 Power supplies range ot (+15 V) — 0— (-9 V) are developed by using thee 
isolated transformers ol rating 230 V / 20-0-15 V (four secondary) and bride rectifiers. 



Fig. 4.3 Application diagram of M57959L 

I'hree lock out circuits are fabricated to get the required blocking period for each 
leg of the inverter. Each lock out circuit consists of two dual mono-stable multivibrators 
IC74123, four AND gate IC 7408, and an HEX inverter IC 7404. The lock out circuit 
provides a time delay of approximately 5 gs between the changes of state in a leg. There 
are 3 possible states in a leg of three-level inverter is shown in table 4.2. The logic 
diagram of lock out circuit is shown in Fig. 4.4. 


Al 

A2 

A3 

A4 

Sal 

Sa2 

1 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

1 

0 


1 

1 

1 

0 


'fable 4.2 Logic for lock out circuit and Encoding Circuit 
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Each phase lequires 4 control signals as there are four switches for each phase. 
An encoder circuit is designed to decrease the no. of output control signals from PC. The 
logic can be obtained by di awing a k-inap from table 4.2. Only two signals are given 
from PC I'oi each phase and they are encoded to four signals. Sal and Sa2 are the two 
signals corresponding to switching signals Al, A2, A3, and A4 given from Digital output 
connector ol PCL 208. The logic diagram for single phase encoding circuit is shown in 
Fig. 4.5. 



Fig. 4.5 Encoding Circuit for one phase 

Encoding circuit for three phases has been designed on single PCB to minimize 
the no. of components. Encoding circuit consists of one NOT gate IC 7404, one EXOR 
gate 1C 7486 and two AND gate ICs 7408. Inputs for encoding circuit are given from 
Digital output iTom PCI. - 208. Output of encoding circuit is given to lock out circuit. 

With the building up of the encoding circuit, lock out circuit and driver circuit, the 
inverter part is ready. I'he dc link voltage is sensed by using a voltage sensor (Model LV 
- 25P). I'wo phase currents are sensed by using two current sensors (Model LA 50P). The 
details ol the voltage and currents are given in Appendix C. An ac tachogenerator is 
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coupled with the induction motor to sense the speed of the motor. A potential divider 
circuit and a suitable filtei aie used to scale down the speed corresponding to suitable dc 
voltage. All the voltage, current and speed sensors are calibrated to determine the 
conveision ol actual and measured quantity. With the building blocks hardware part of 
the implementation is over. 

4.3 PC - based Implementation of the Control Circuit 

I he control circuit mainly consists of a PC operating at a clock frequency of 266 
MHz. I he control algorithm executed in a real time on the PC. The calculations are 
performed in.sidc the PC. A high speed Data Acquisition Card (PCL - 208) has been 
interfaced with the PC for the data transfer between PC and external hardware. The 
detailed specification of the card is given in Appendix B. In the actual experimental 
setup, the drive is controlled by a software program. In the PC-based control scheme the 
data, like motor currents, motor speed, and dc link voltage signals are given to the digital 
computer through the A/D converter of high speed data acquisition card (PCL - 208). 
The switcliing status for inverter are generated by calculations performed inside the 
computer and control pulses, for inverter gate circuits are computed by the digital output 
ports of the data acquisition card. 

4.4 PC - based Control of the Proposed DTFC Scheme 

The PC-based real time program is given in Appendix D. the program is in C 
language. In the control program, analog signals like motor phase currents, motor speed, 
and dc link voltage are sensed. In fact the sensed analog signals are converted into digital 
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signals by the in-built A/D converter of the DAC card. In the algorithm the converted 
digital values will again be converted to analog values. The converted analog values after 
multiplication with sensor gains will give the actual values of the voltage, currents, and 
speed piesent at that time. These actual analog values will be used by the DTFC 
algoiithm to calculate the sector number and the status of the flux and torque. From the 
knowledge of these tin ee variables, namely flux status, torque status and sector angle the 
required voltage vector and hence switching signals for inverter are selected from the 
switching table. Finally the PC will send out six switching signals Sal, Sa2; Sbl, Sb2; 
and Scl, Sc2; two tor each phase, required to operate the machine in direct torque and 
flux control strategy. 

4.5 Expcriinentiil Results 

The experiment is done on a 3-pahse, 1.1 kW, 415 V, 2.6 A, 1410 rpm and 50 Hz, 
star connected squirrel cage induction motor which is mechanically coupled with a 230 
V, 1500 rpm dc machine for loading purpose. Typical data of machine parameters are 
given in Appendix ‘A’. Typical experimental results are presented and they are compared 
with the simulation results to validate the design methodology. 

4.5.1 Under no load condition 

Fig. 4,6 - 4.8 show the simulation and experimental results under no load 
condition. F'ig 4.6 shows that for the step change of speed the stator flux is constant and is 
equal to the reference value of 1.07 Wb. Fig. 4.7 shows that the torque developed when 
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there is a step change in speed in the starting. Fig. 4.8 shows the cuirent waveform under 

step change in speed. Fig. 4.7 and Fig. 4.8 have been compared with simulation results. 

4.5.2 Speed reversal 

Simulaltd cind txptnmental waveforms of the speed response, the flux response, 
motor phase current and torque response for the step change in speed from 280 ipm 
(forward motoring) to -280 rpm (reverse motoring) and vise versa are shown in Fig. 4.9 
to 4.1 1. The speed response is found to be fast. During the transient period, the motor 
phase current is high as shown in Hg. 4.1 1. The flux is constant during the transient and 
is equal to the reference value 1.07 Wb. The torque response is shown in the Fig. 4.10 
during speed reversal. This shows the de-coupling of the torque and flux and establishes 
that the operation i.s .similar to that of separately excited dc machine. 

4.5.3 Step change in torque command 

'fhe load on induction motor is applied by connecting a resistive load to the dc 
generator, which is coupled to the motor. To bring about a step change in torque 
command, dc generator coupled to the induction motor was loaded by switching on the 
field supply while the resistive load is connected to the dc generator. In Fig. 4.12, it is 
shown that as the load is applied on motor, flux is constant and is at reference value. Fig. 
4.13 show.s that llie speed is almost constant at its reference value when there is an 
application of load torque. Fig. 4.14 shows the current waveform under step change in 
load torque. 
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4.6 Conclusion 

Expel iniental setup for the direct torque and flux control of induction motor is 
explained in this chaptei . Data transfer between PC and hardware is achieved by PCL — 
208 data acquisition card. Control signals are generated in side the computer. Instead of 
conventional designed control circuits, a high performance M59759L hybrid drives ICs 
are used tor amplification of control pulses and to provide isolation between low power 
control and high power control and high power IGBT inverter circuit, thus making 
inverter design compact. The simulation and experimental results are compared in this 
chapter and experimental results are matched with simulated results. The experimental 
results for bolli steady state and transient conditions were found satisfactory. The torque 
response is almost instantaneous and speed of the motor is fast. From the result it is seen 
that the machine is under decoupled control, where torque and flux can be controlled 
directly and independently. 
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1 div = 0.65wb (Stator flux) 
Fig. 4.6(a) Simulated Waveform 


Tek JU •Stop MPosi-smoim 

.r'* 


SASfEmic 


setups 


SMJfce 



P.ef 

i 


Sdve 


flif 

gi 

CHI 1.0CIV CH2 50.0mV M 1.00s CHI y O.OCV 

Ref A 1.00 V 1.00s RefB 50.0mV I.Ws < 10 H 2 


X-axis: Idiv = 1 sec (time) 

Y-axis: Idiv = 2Wb (stator flux, Channel 1) 

: Idiv = 230 rpm (Motor Speed, Channel 2) 

(There is a, +2.5V offset in D/A converter’s output, Channel 1) 

Fig. 4.6 Flux response during step change of speed from 0 rpm to 280 rpm 
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1 div = 2 N-m (Torque developed) 


Fig. 4.7(a) Simulated Waveform 
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Fig. 4.7 Torque response during step change in speed from 0 rpm to 280 rpm 
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1 div = 4.5 A (Motor phase ‘a’ current) 


Fig. 4.8 (a) Simulated Waveform 
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Fig. 4.8 Motor phase ‘a’ cuixent during step change in speed from 0 rpm to 280 rpm 
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1 div = 2wb (Stator flux) 


Fig. 4.9 (a) Simulated Waveform 
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Fig. 4.9 Stator flux response when step change in speed from 280 rpm (forward 
motoring) to -280 rpm (reverse motoring) and vise versa. 
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1 div = 7.25 N-m (Torque developed) 
Fig. 4.10(a) Simulated Waveform 
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Fig. 4.10 Torque response when step change in speed from 280 rpm (forward motoring) 
to -280 rpm (reverse motoring) and vise versa. 
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1 div = 4.5 A (Motor phase ‘a’ current) 


Fig. 4.1 1(a) Simulated Waveform 
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Fig. 4.1 1 Motor phase ‘a’ current when step change in speed from 280 rpm (forward 
motoring) to -280 rpm (reverse motoring) and vise versa. • 
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Fig. 4.12(a) Simulated Waveform 
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Fig. 4.12 Stator flux response tmder step change in load 
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1 div = 200 rpm (Motor Speed) 


Fig. 4.13(a) Simulated Waveform 
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Fig. 4.13 Speed response under step change in load 


Chapter 4: PC — based Implementation of the Induction motor drive 


84 




1 div = 3 A (Motor phase ‘a’ current) 
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Fig. 4.14(a) Simulated Waveform 
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Fig. 4.14 Motor phase ‘a’ current under step change in load 
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X-axis 1 div = 0.2 sec (time) 
Y-axis 1 div = O.Swb (Stator flux) 


Fig. 4.15(a) Simulated Waveform 
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Fig. 4.15 d and q-axis fluxes, under steady state running condition. 
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A View of Three-Level Inverter with Lock-out Circuit 
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A Inside View of Lockout Circuit 
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A View of DTFC control of Induction Motor Drive System under operation 



A Close View of Sensor Circuits and Encoder Circuit 
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Chapter 5 


Conclusions 


5.1 Contributions of the Present Thesis work 

Initially, simulation of four-quadrant direct torque and flux control of induction 
motor has been presented. The results show that the harmonics in the source current is 
less, and power factor is unity. Subsequently, PC-based implementation of a direct torque 
and flux controlled three-level inverter fed induction motor has been presented for high 
power applications. The dynamic response of the drive is fast with direct and independent 
control of torque and flux. From the simulation results it can be seen that the low speed 
performance of the induction motor has been improved using modified look up table. 

The salient contributions of the present thesis work can, thus, be summarized as 
follows. 

1 . Simulation of three-level synchronous link converter. 

2. Modeling and simulation of a utility friendly direct torque and flux controlled 
induction motor drive system. 

3. PC-based hardware implementation of the proposed, direct torque and flux 
control scheme, for three-level inverter-fed induction motor. 



5.2 Scope for the Future Work 

In the present work hysteresis current control scheme is used for SLC. For high 
power application, hysteresis current control scheme cannot be used. This is due to the 
switching limitation of high power devices. Constant frequency pulse width modulated 
control scheme is suitable for high power application. 

Direct torque control (DTFC) of induction machines has been actively researched 
in the area of AC machine drives for a number of years. DTFC method described in this 
work relies on the hysteresis principle in which both torque values and stator flux 
magnitudes are regulated within pre-specified numerical ranges. The performance of such 
a scheme depends on the error band set between the desired and measured torque values. 
In general the narrower this band, the superior the torque performance, however the 
higher the switching frequency. The flux and torque hysteresis bands are the only gains to 
be adjusted in DTFC, and the switching frequency and the current waveform are greatly 
influenced by them. Therefore, the magnitude of the hysteresis band should be 
determined based on reasonable guidelines which can avoid excessive switching 
frequency and current harmonics in the whole operating region. 

In the present thesis, the reference torque is generated by using the motor speed 
signal obtained from a mechanical transducer mounted on the rotor shaft. However, this 
transducer makes the system less robust and reduces the system reliability. An improved 
drive system can be designed by estimating the speed to avoid the use of speed sensor. 

In the proposed scheme, change in parameters is not taken into consideration. If 
the parameter variation is also considered the accuracy of the results can be improved. 
Thus if an adaptive motor model is taken the problem of change in characteristics with a 
change in machine parameters can be tackled and results can be improved. 
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Two major problems that are usually associated with DTFC drives are: 1) 
switching frequency that varies with operating conditions and 2) high torque ripple. To 
solve these problems, a constant switching frequency torque controller has to be used 
instead of the conventional hysteresis-based controller. 
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Appendix A 

Specifications of a 3 - phase Induction Motor 


Ratings 


Power ; 

1.1 kW 

Voltage : 

415 L-L 

Current : 

2.6 A 

Speed : 

1410 rpm 

Moment of Inertia of motor-generator set (j) : 

0.012 kg-m' 

Stator Parameters 

Resistance (rs) : 

8.2 Ohm 

Leakage inductance (11s) : 

0.026 H 

Magnetizing inductance (Im) : 

0.377 H 

Rotor Parameters (ref. to stationary) 

Resistance (rr) : 

: 4.3 Ohm 

Leakage inductance (Hr) 

: 0.026 H 

Specifications of DC machine 

Ratings 

Power 

: 1 hp 

Voltage 

; 230 V 

Current 

: 4.5 A 

Speed 

: 1500 rpm 



Appendix B 


Specifications of PCL — 208 Data Acquisition Card 

PCL - 208 is a high performance, high speed multifunction data acquisition card 
for the IBM PC/XT/AT or compatibles. The high end specifications of this full size card 
and complete software support make it deal for wide range of applications in the 
industrial and laboratory environment, like data acquisition, process control, automatic 
testing and factory automation. 

Main Features: 

> Switch selectable 1 6 single-ended or 8 analog input channels. 

> An industrial standard 12 bit successive approximation converter (ADC 674) to 
convert analog inputs. The maximum A/D sampling rate is 60 kHz in DMA 
mode. 

> Switch selectable versatile analog input ranges. 

Bipolar: +/- 0.5 V, +/- 1 V, + 1 - 1.5 V, +/- 5 V, +/- 10 V. 

Unipolar: +1 V,+2 V,+5 V,+10 V. 

> Provides three A/D trigger modes: software trigger. Programmable pacer trigger 
and external trigger pulse trigger. 

> A/D converted data can be transferred by program control, interrupt handler 
routine or DMA transfer. 

> An INTEL 8254 programmable Timer/Counter provides pacer output (trigger 
pulse) at the rate of 2.5 MHz to 71 minutes/pulse to the A/D. The timer time base 
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is switch selectable to 10 MHz or 1 MHz. One 16-bit counter channel is reserved 
for user configuration applications. 

y Two 12-bit monolithic multiplying D/A output channels. Output range of 0 to +5 
V can be created by using the onboard -5 V reference. This precision reference is 
derived from the A/D converter reference. External AC or DC reference can also 
be used to generate other D/A output ranges. 

> TTL/DTL compatible 16 digital input & 16 digital output channels. 

A/D & D/A Converter Specifications 
Analog Input Specifications 

o Channels : 16 Single-ended or 8 Differential switch selectable, 

o Resolution : 12 bits. 

o Input range : Unipolar: +1 V,+2 V,+5 V,+10 V. 

Bipolar: +/- 0.5V, +/- IV, +/-2.5V, +/- 5V, +/- lOV. 
All input ranges are switch selectable, 
o Over voltage : Continuous +/- 30 V Max. 

o Conversion Type : Successive Approximation, 

o Conversion Speed : 60 kHz max. 

o Accuracy Speed : 0.01 % of reading +/- 1 bit. 

o Linearity : +/- 1 bit. 

o Trigger Mode : Software trigger, onboard programmable timer or 

external trigger. 

o Data transfer : Program control. Interrupt control or DMA. 
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Analog Output Specifications: 


o Channels : 

o Resolution : 

o Output range : 

o Reference voltage : 

o Conversion Type : 

o Linearity : 

o Output Drive : 

o Settling Time : 

General Specifications: 
o Power Consumption : 


o I/O Connector 


o I/O Base Address 


2 channels. 

12 bits. 

0 to +5 V with fixed -5 V refemce. 

+/- 10 V with external DC to AC reference. 

Internal: -5 V (+/- 0.05 V). 

External: DC or AC, +/- 10 V max. 

12 bit monolithic multiplying (DAC 7541) 

+/- 1/2 bit. 

+/- 5ma max. 

5 microseconds. 

+5 V : typ.700 mA, max. lA. 

+12 V : typ.140 mA, max. 200 mA. 

-12 V : typ.l4 mA, max. 20 mA. 

20 pin flat cable connector for all Analog/Digital 
I/O ports. 

Requires 16 consecutive address locations. Base 
address is definable by the DIP switches for address 
lines A9 - A4. 


Specifications used in Experimental Work: 

• 16 Single-ended Analog input channels 

• Bipolar +/- 1 0 V, -5 V internal reference. 
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Appendix C 

Voltage and Current Sensors 

Technical Specifications of Voltage Sensors, Current Sensors, are tabulated in the 
following tables. 


Table C. 1 : Voltage Sensor: Model LV 25 - P 


Primary nominal current 

: 10 mA 

Secondary nominal current 

: 25 mA 

Conversion ratio 

: 2500:1000 

Overall accuracy at 25° C 

: ±0.9 % 

Supply Voltage 

: ±12 V (±5%) 

Isolation test voltage 

: 2.5 kV (rms) / 50 Hz / 1 min 

Linearity 

: <0.2 % 

Response time 

: <40 ps 

Operating Temperature 

: 0° C to 70° C 

Current consumption 

: 1 0 mA ± output current 

Primary Internal Resistance 

: 250 Q (at 70° C) 

Secondary Internal Resistance 

: 110 0(at70°C) 

Weight 

: 22 gm 

Operating Range 

: 10 to 500 V 


Polarity Marking: A positive output current is obtained on terminal M when positive 
voltage is applied in terminal +HT of the primary circuit. 
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Primary resistor R1 : The transducer’s optimum accuracy is obtained with the nominal 
primary current. As far as possible, Ri should be calculated so that the nominal voltage to 
be measured corresponding to primary current of 10 mA. 


Table C.2: Measuring Resistances 


Voltage rating 

Rm (min) 

Rm (max) 

With ±12 V and ±10 niA max 

30 n 

190 a 

With ±12 V and ±14 mA max 

30 n 

100 


Table C.3: Connection Pins 


Pin ± 

: Supply V oltage ± 1 5 V 

Pin M 

: Measuring point 

Pin - 

; Supply Voltage - 1 5 V 

Pin ±HT 

: Primary Voltage ± 

Pin-HT 

: Primary Voltage - 
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Table C.4; Current Sensor: Model LA 55 - P 


Primary nominal cunent 

: 50 A 

Secondary nominal current 

: 50 mA 

Conversion ratio 

: 1:1000 

Overall accuracy at 25® C 

: ±0.65 % 

Supply Voltage 

: ±12 V (± 5%) 

Isolation test voltage 

: 2.5kV(rms)/50Hz/ 1 min 

Linearity 

: <0.15% 

Response time 

: <1 |is 

Operating Temperature 

: -25° C to 85° C 

Current consumption 

: 1 0 mA + output current 

Weight 

: 18 gm 


Table C.5: Measuring Resistances (85® C) 


Voltage rating 

Rm (min) Rm (max) 

With ±12V and ±50mA max 

60 0 

95 0 

With ±12V and ±70mA max 

60 0 

60 0 

Table C.6: Connection Pins 


Pin ± : 

Supply Voltage +15 V 


PinM : 

Measuring point 


Pin- : 

Supply Voltage -15 V 
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Appendix D 


1. Simulation program 

#include<stdio .h> 

#include<conio .h> 

#include<stdlib.h> 

#include<math.h> 

#defme pi 3.1415 
void mainQ 
{ 

int p=4,f=50,x=100; 

float vrated=41 5.0, vqs=0.0,vds=0.0, pm; 

float wrm; 

float h=0.00001; 

float Im=0.57,lr=0.632,ls=0.632,rs=l .0,rr=3.46; 
float 11s, Hr, 1ml; 
float j=0.02,b=0.01; 

float ids=0.0,iqs=0.0,idr=0.0,iqr=0.0,is=0.0,wr=0.0,wrl=0.0; 

float fds=0.0,fqs=0.0,fdr=0.0,fqr=0.0,fs=0.0; 

float fds 1 =0. 0,fqs 1 =0.0,fdr 1 =0.0,fqr 1 =0.0; 

float wre=0.0,wre0=0.0; 

float wr_ref, wrm_ref=0 .0 ,wr_rated,wc; 

float fqm=0.0,fdm=0.0; 

float te=0.0,tmx,tmn,ted,te_ref=0.0,tl=0.0,tl_rated=24.536; 
float fref,fh,fl; 

float tetae=0 . 0,ki= . 08 ,kp= . 3 222 ,ht 1 ,ht2,hf; 

float sim_time=5.0,t,ti,k; 
int i,ct=0,cf=0,tetan=0,v; 

float kl [4],k2[4],k3[4],k4[4],k5[4]; 

// Variables for SLC 

float vdc_ref=280.0,vdc=0.0,vm=280.0,vas=0.0,vbs=0.0,vcs=0.0; 

float vdcl=0.0,vdc2=0.0,nr=1410.0; 

float van=0 . 0 , vbn=0 . 0 , vcn=0 . 0,vao=0 .0 ,vbo=0 . 0,vco=0 . 0 ; 

float las_ref=0.0,lbs_ref=0.0,lcs_ref=0.0,las=0.0,lbs=0.0,lcs=0.0; 

float Im=0.0,Idcl=0.0,Idc2=0.0; 

float lsr=.00024,Rr=35.5,c=0.008,kpr=.2,kir=l . 1 95,il=0.0; 

float v_err=0.0,err0=0.0,errl=0.0,sal=0.0,sa2=0.0,sbl=0.0,sb2=0.0,scl=0.0,sc2=0.0; 
float wrs=2.0*pi*50.0; 
float hi=0.2; 

float kcl [4],kc2[4],kc3[4],kc4[4],kc5[4]; 



float Ial=0.0,Ibl=0.0,Icl=0.0,vdl=0.0,vd2=0.0; 

/* Lookup Table */ 
const int lookuptable[2][5][12]={ 

{ 

{52,53,62,63,12,13,22,23,32,33,42,43}, 

{51,53,61,63,11,13,21,23,31,33,41,43}, 

{8,7,8,7,8,7,8,7,8,7,8,7}, 

{31,33,41,43,51,53,61,63,11,13,21,23}, 

{32,33,42,43,52,53,62,63,12,13,22,23}, 

{ 

{62,63,12,13,22,23,32,33,42,43,52,53}, 

{61,63,11,13,21,23,31,33,41,43,51,53}, 

{7,8,7,8,7,8,7,8,7,8,7,8}, 

{21,23,31,33,41,43,51,53,61,63,11,13}, 

{22,23,32,33,42,43,52,53,62,63,12,13} 

} 

}; 


FILE *fpl,*fp2,*fp3,*fp4,*fp5,*fp6,*fp7,*fp8,*fpl0,*fpll,*fpl2,*fpl3,*fpl4,*fpl5,*fpl6; 

fpl=fopen("3torql.txt","w"); 
fp2=fopen("3spdl.txt","w"); 
fp3=fopen("3Iqsl.txt","w"); 
fp4=fopen("3fqsl .txt","w"); 
fp5=fopen("3fluxl .txt","w"); 
fp6=fopen("3fdsl.txt","w"); 
fp7=fopen("3fqdsl.txt","w"); 
fp8=fopen("3Idsl.txt","w"); 

//fp9=fopeii("3Comd.txt","w"); 
fp 1 0=fopen("3testl .txt","w"); 
fp 1 1 =fopen(" 3 Vlas 1 .txt" , " w"); 
fpl2=fopen("3VIbsLtxt","w"); 
fpl 3=fopen("3 Vies 1 .txt'V'w"); 
fp 1 4=fopen("3 Vdc 1 .txt","w"); 
fpl 5=fopen("3Iabcsl. txt'V'w"); 
fpl6=fopen("3Vdcl2Ltxt","w"); 

clrscrO; 

printf("Program is running... "); 


lls=ls-lm; 

llr=lr-lni; 
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Im 1 = 1 /(( 1 /lm)+( 1 /lls)4-( 1 /Hr)); 


fref=sqrt(2 . 0/3 . 0) * (vrated/wrs) ; 

htl=0.01*tl_rated; 

ht2=0 . 02 * tl_rated; 

hf=0.01*fref; 

wr_rated=nr*2.0*pi/60.0; 

tmx=2.0*tl_rated; 
tmn=-2 . 0 * tl_rated ; 

for (t=0.0;t<=sim_time;t+=h) 

{ 

if(t<0.5) 

wr_ref=0.05 * wr_rated; 
if(t>=.50) 

wr_ref=0.53 19*wr_rated; 
if((t>=1.75)&&(t<3.0)) 
wr_ref=-0.53 1 9*wr_rated; 

//tl=5.0; 

if(t>=3.0) 

wr_ref=0 .531 9 *wr_rated; 

//tl=0.0; 

wre=wr_ref-wr; 

te_ref=te_ref+kp*(wre-wreO)+ki*h*wre; 

if(te_ref>=tmx) 

te_ref=tmx; 

if(te_ref<=tmn) 

te_ref=tmn; 

wre0=wre; 

ted=te_ref-te; 

if(ted>=0.0 && ted<htl && ct==0) 
ct=0; 

else if(ted>=htl && ted <ht2 && ct==0) 
ct=l; 

else if(ted>=ht2) 
ct=2; 

else if(ted>=htl && ted<ht2 && ct==2) 
ct=2; 

else if(ted<htl && ted>=0.0 &&ct==2) 
ct=l; 

else if(ted<=0.0 && ct==l) 
ct=0; 

else if(ted<=0.0 && ted>-htl && ct=0) 



ct=0; 

else if(ted <= -htl && ted>-ht2 && ct==0) 
ct=-l; 

else if(ted<=-ht2) 
ct=-2; 

else if(ted<= -htl && ted>-ht2 && ct=-2) 
ct=-2; 

else if(ted>-litl && ted<=0.0 &&ct==-2) 
ct=-l; 

else if(ted<0.0 && ted>=-htl && ct==-l) 
ct=-l; 

else if(ted>=0.0 && ct==-l) 
ct=0; 

// else printf("ct=%d\n",ct); 

fh=fref+hf; 

fl=fref-hf; 

if(fs<=fl) 

cf=l; 

if(fs>=fh) 

cf=0; 


// Calculation of flux vector 

if(tetae>=0.0 && tetae<pi/6.0) 
tetan=l; 

else if(tetae>=pi/6.0 && tetae<pi/3.0) 
tetan=2; 

else if(tetae>=pi/3.0 && tetae<pi/2.0) 
tetan=3; 

else if(tetae>=pi/2.0 && tetae<2.0*pi/3.0) 
tetan=4; 

else if(tetae>=2.0*pi/3.0 && tetae<5.0*pi/6.0) 
tetan=5; 

else if(tetae>=5.0*pi/6.0 && tetae<pi) 
tetan=6; 

else if(tetae>=pi && tetae<7.0*pi/6.0) 
tetan=7; 

else if(tetae>=7.0*pi/6.0 && tetae<4.0’''pi/3.0) 
tetan=8; 

else if(tetae>=4.0*pi/3.0 && tetae<3.0*pi/2.0) 
tetan=9; 

else if(tetae>=3.0*pi/2.0 && tetae<10.0*pi/6.0) 
tetan=10; 

else if(tetae>=10.0*pi/6.0 && tetae<11.0*pi/6.0) 
tetan=ll; 



else tetan=12; 


// Syncronous link converter 
err 1 =vdc_ref- vdc ; 

v_err=v_err+kpr*(errl -errO)+kir*h*errl ; 

errO=errl; 

if(v_err>5.0) 

Im=5.0; 

else if(v_err<-5.0) 

Im=-5.0; 
else Im=v_err; 

1 as_re 1 = I m * s i n (wrs * t) ; 

Ibs_ref=Im * sin(wrs * t-2 .0 *pi/3 . 0) ; 
Ics_ref=Im*sin(wrs*t-4.0*pi/3.0); 

if(Ias_ref-Ias<-hi) {sal=l .0;sa2=0.0;} 
else if(Ias_ref-Ias>hi) {sal=0.0;sa2=1.0;} 
else {sal=0.0;sa2=0.0;} 

if(Ibs_ref-Ibs<-hi) {sbl=1.0;sb2=0.0;} 
else if(Ibs_ref-Ibs>hi) {sbl=0.0;sb2=1.0;} 
else {sbl=0.0;sb2=0.0;} 

if(Ics_ref-Ics<-hi) {scl=l .0;sc2=0.0;} 
else if(Ics_ref-Ics>hi) {scl=0.0;sc2=1.0;} 
else {scl=0.0;sc2=0.0;} 

vas=vm* sin(wrs* t) ; 

vbs=vm*sin(wrs*t-2.0*pi/3.0); 

vcs=vm*sin(wrs*t-4.0*pi/3.0); 

vao=sa 1 * vdc 1 -sa2 * vdc2 ; 
vbo=sb 1 * vdc 1 -sb2 * vdc2 ; 
vco=scl *vdcl-sc2*vdc2; 

van=( 1. 0/3.0)* (2.0* vao- vbo- vco) ; 
vbn=(1.0/3.0)*(2.0*vbo-vao-vco); 
vcn=(1.0/3.0)*(2.0*vco-vao-vbo); 

for(i=0;i<4;++i) 

{ 

kc 1 [i]=h*(vas-van-Rr*Ial )/lsr; 
kc2 [i] =h* (vbs- vbn-Rr * Ib 1 )/l sr ; 



kc3 [i]=h* (vcs-vcn-Rr* Ic 1 )/lsr ; 

kc4 [i]=h* (sa 1 * la 1 +sb l*Ibl+scl*Icl -il)/c; 

kc5 [i]=-h* (sa2*IaH-sb2*Ib 1 +sc2*Ic 1 +il)/c; 

if(i==Oi|i==l) 

{ 

Ial+=kcl[i]/2.0; 

Ibl+=kc2[i]/2.0; 

Icl+=kc3[i]/2.0; 

vdl+=kc4[i]/2.0; 

vd2+=kc5[i]/2.0; 

} 

if(i==2) 

{ 

Ial+=kcl[i]; 

IbH-=kc2[i]; 

Icl+=kc3[i]; 

vdl+=kc4[i]; 

vd2+=kc5[i]; 

} 

} 

Ias+=(kcl[0]+2.0*kcl[l]+2.0*kcl[2]+kcl[3])/6.0; 

Ibs+=(kc2[0]+2.0*kc2[l]+2.0*kc2[2]+kc2[3])/6.0; 

Ics+=(kc3[0]+2.0*kc3[l]+2.0*kc3[2]+kc3[3])/6.0; 

vdcl+=(kc4[0]+2.0*kc4[l]+2.0*kc4[2]+kc4[3])/6.0; 

vdc2+=(kc5[0]+2.0*kc5[l]+2.0*kc5[2]+kc5[3])/6.0; 

Ial=Ias; 

Ibl=Ibs; 

Icl=Ics; 

vdl=vdcl; 

vd2=vdc2; 

Idc 1 =sal * las+sb 1 * Ibs+sc 1 * Ics; 

Idc2=sa2 * Ias+sb2 *Ibs+sc2*Ics; 
vdc=vdcl+vdc2; 

//synchronous link converter end 

// Selection of voltage vector from look up table 

v=lookuptable[cf] [ct+2] [tetan- 1 ] ; 
switch(v) 

{ 

case 1 1 : 

vds=(1.0/3.0)*vdc; 



vqs=0.0; 
break; 
case 12: 

vds=(2.0/3.0)*vdc; 
vqs=0.0; 
break; 
case 2 1 : 

vds=(1.0/6.0)*vdc; 
vqs=vdc/(2.0* 1.732); 
break; 
case 22: 

vds=(1.0/3.0)*vdc; 
vqs=vdc/( 1.732); 
break; 
case 3 1 : 

vds=-(l. 0/6.0)* vdc; 
vqs=vdc/(2.0* 1.732); 
break; 
case 32: 

vds=-(l-0/3.0)*vdc; 
vqs=vdc/1.732; 
break; 
case 41: 

vds=-(1.0/3.0)*vdc; 
vqs=0.0; 
break; 
case 42: 

vds=-(2.0/3.0)*vdc; 
vqs=0.0; 
break; 
case 5 1 : 

vds=-(1.0/6.0)*vdc; 
vqs=-vdc/(2.0* 1 .732); 
break; 
case 52: 

vds=-(1.0/3.0)*vdc; 
vqs=-vdc/(1.732); 
break; 
case 61: 

vds=(1.0/6.0)*vdc; 
vqs=-vdc/(2.0* 1.732); 
break; 
case 62: 

vds=(1.0/3.0)*vdc; 

vqs=-vdc/1.732; 

break; 
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case 13: 

vds=(1.0/2.0)*vdc; 
vqs=vdc/(2.0* 1.732); 
break; 
case 23 : 

vds=0.0; 
vqs=vdc/(1.732); 
break; 
case 33: 

vds=-(1.0/2.0)*vdc; 
vqs=vdc/(2.0* 1.732); 
break; 
case 43 : 

vds=-(1.0/2.0)*vdc; 
vqs=-vdc/(2.0* 1.732); 
break; 
case 53: 

vds=0.0; 

vqs=-vdc/(1.732); 
break; 
case 63 : 

vds=(1.0/2.0)*vdc; 
vqs=-vdc/(2.0* 1.732); 
break; 
case 7,8: 
default: 

vds=0.0; 

vqs=0.0; 

break; 

} 

// printf("%fyt%f\n",vds,vqs); 

// fprintf(fp9,"%0.4f\t%0.5f\t%d\t%d\t%d\t%d\t%0. lf\t%0. 1 f\n",t,tetae,cf,ct,tetan,v,vds,vqs); 
for(i=0;i<4;++i) 

{ 

kl [i]=h*(vqs-(rs/lls)*(fqsl-((lml/lls)*fqsl + (lml/llr)*fqrl))); 
k2[i]=h*(vds-(rs/lls)*(fdsl-((lml/lls)*fdsl + (lml/llr)*fdrl))); 
k3[i]=h*(wrl*fdrl-(rr/llr)*(fqrl-((lml/lls)*fqsl + (linl/llr)*fqrl))) ; 
k4[i]=h*(-wrl*fqrl-(rr/llr)*(fdrl-((lml/lls)*fdsl + (lml/llr)*fdrl))) ; 
k5[i]=h*(p/(2.0*j))*(te-tl-b*(2.0/p)*wrl); 


if(i==0||i==l) 

{ 

fqsl+=kl[i]/2.0; 

fdsl+=k2[i]/2.0; 

fqrl+=k3[i]/2.0; 

fdrl+=k4[i]/2.0; 


no 



wrl+=k5[i]/2.0; 

} 

if(i==2) 

{ 

fqsl+=kl[i]; 

fdsl+=k2[i]; 

fqrl+=k3[i]; 

fdrl+=k4[i]; 

wrl+=k5[i]; 

} 

} 

fqs+=(kl [0]+2.0*kl [l]+2.0*kl [2]+kl [3])/6.0; 
fds+=(k2[0]+2.0*k2[l]+2.0*k2[2]+k2[3])/6.0; 
fqr+=(k3 [0]+2.0*k3 [1 ]+2.0*k3 [2]+k3 [3])/6.0; 
fdr+=(k4[0]+2.0*k4[l]+2.0*k4[2]+k4[3])/6.0; 
wr+=(k5[0]+2.0*k5[l]+2.0*k5[2]+k5[3])/6.0; 

fs=sqrt(fds*fds+fqs*fqs); 

fqsl=fqs; 

fdsl=fds; 

fqrl=fqr; 

fdrl=fdr; 

wrl=wr; 

fqm=((lml/lls)*fqs + (lml/llr)*fqr); 
fdm=((lml/lls)*fds + Oml/llr)*fdr); 

iqs=(fqs-fqm)/lls; 

ids=(fds-fdm)/lls; 

iqr=(fqr-fqm)/llr; 

idr=(fdr-fdm)/llr; 

te=(3.0/2.0)*(p/2.0)*(Ms*iqs-fqs*ids); 
is=sqrt(ids*ids+iqs*iqs); 
wrm=wr*2.0/p; 
wrm_ref=( wr_ref‘ 2 . 0/p) ; 

pm=te*wrm+3.0*is*is/2.0*rs; 

il=pin/vdc; 

if(fqs==0 . 0&&fds==0 . 0) 
tetae=0.0; 
else 

tetae=atan2(fqs,fds); 
if(tetae<0.0) 
tetae= tetae+2.0*pi; 



if(x>=100) 

{ 

fprintf(fp 1 , " %f\t%f\t%f\n" 

fprintf(fp2,"%f\t%f\t%0.4f\n",t,wr*60.0/(2.*pi),wr_ref‘60.0/(2.*pi)); 

fprintf(fp3,"%f\t%f\n",t,iqs); 

fprintf(fp4,"%f\t%f\n",t,fqs); 

fprintf(fp5,"%f\t%f\t%f\n",t,fs,fref); 

fprintf(fp6,"%f\t%f\n",t,fds); 

fprintf(fp7,"%f\t%f\n",fqs,fds); 

fprintf(fp8,"%f\t%An",t,ids); 

fprintf(fplO,"%f\t%f\t%f\t%f\t%f\t%f\t%f\n",t,il,pm,Im,v_err,Ias,Ias_ref); 

fprintf(fp 1 1 ,"%f\t%f\t%f\n",t,vas/l 0.0, las); 

fprintf(fpl 2,"%f\t%f\t%f\n",t,vbs/l 0.0, lbs); 

fprintf(fpl 3,"%f\t%f\t%f\n",t,vcs/l 0.0,lcs); 

fprintf(fpl4,"%f\t%f\t%f\n",t,vdc_ref,vdc); 

fprintf(fpl5,"%f\t%f\t%f\t%f\n",t,Ias,Ibs,Ics); 

fprintf(fpl 6,"%f\t%f\t%f\n",t,vdcl ,vdc2); 

x=0; 

} 

X++; 

} 

fclose(fpl); 

fclose(fp2); 

fclose(fp3); 

fclose(fp4); 

fclose(fp5); 

fclose(fp6); 

fclose(fp7); 

fclose(fp8); 

// fclose(fp9); 
fclose(fpll); 
fclose(fpl2); 
fclose(fpl3); 
fclose(fpl4); 
fclose(fpl5); 
fclose(fpl6); 


} 



2. Listing of Real Time PC-Based Program 

#include<stdio.h> 

#include<conio .h> 

#include<stdlib.h> 

#include<math.h> 

#include <dos.h> 

#include<time.h> 

#define pi 3.14159 
#define BASE 0x300 
#define TRUE 1 
void main() 

{ 

int p=4,f=50; 

float vdc,vrated=4 1 5 .0,vqs=0.0,vds=0.0; 
float wrs,wrm; 
float h=0.0002; 

float lm=0.377,lr=0.403,ls=0.403,rs=1.0,rr=4.3; 
float 11s, Ur, 1ml; 
float j=0.0 1 2,b=0.0; 

float ids=0.0,iqs=0.0,idr=0.0,iqr=0.0,is=0.0,wr=0.0,wrl=0.0; 
float fds=0.0,fqs=0.0,fdr=0.0,fqr=0.0,fs=0.0; 
float fdsl=0.0,fqsl=0.0,fdrl=0.0,fqrl=0.0; 
float wre=0.0,wre0=0.0; 

float wr_ref,wrm_ref=0.0,wr_rated,nr=1440.0,wc; 
float fqm=0.0,fdm=0.0; 

float te=0.0,tmx,tmn,ted,te_ref=0.0,tl=0.0,tl_rated=24.536; 
float fref,fh,fl; 

float tetae=0.0,ki=l .08,kp=2.5,htl,ht2,hf; 
float t,k; 

int i,ct=0,cf=0,tetan=0,v; 

float Id [4],k2[4],k3[4],k4[4],k5[4]; 

int muxscan=0x30,swlogic; 
int trigger=0x00; 
int check, u; 

int mb,lb,hb,msb[4],lsb[4]; 
float wra,va=0.0,vad=0.0,la,lb,lc; 
float val[4]; 

clock_t strt=0,end=0; 

int sal=l,sa2=l,sbl=l,sb2=l,scl=l,sc2=l; 
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/* Lookup Table */ 
const int lookuptable[2][5][12]={ 

{ 

{ 52,53 ,62,63, 1 2,13 ,22,23,32,33,42,43 }, 
{51,53,61,63,11,13,21,23,31,33,41,43}, 
{8,7,8,7,8,7,8,7,8,7,8,7}, 

(3 1 ,33,41 ,43,5 1 ,53 ,6 1 ,63,1 1 ,1 3,2 1 ,23 }, 
{ 32,3 3 ,42,43 ,52,53 ,62,63 , 1 2, 1 3 ,22,23 } , 

}, 

{ 

{62,63,12,13,22,23,32,33,42,43,52,53}, 

{61,63,11,13,21,23,31,33,41,43,51,53}, 

{7,8,7,8,7,8,7,8,7,8,7,8}, 

{21,23,31,33,41,43,51,53,61,63,11,13}, 

{22,23,32,33,42,43,52,53,62,63,12,13} 

} 

}; 


wrs=2.0*pi*f; 

lls=ls-lm; 

llr=lr-lm; 

Im 1 = 1 /(( 1 /lm)+( 1 /lls)+( 1 /Ur)); 

fref=sqrt(2 .0/3 . 0) * (vrated/wrs) ; 
htl =0.00 1 *tl_rated; 
ht2=0.002*tl_rated; 
hf=0.01*fref; 
wr_rated=nr*2.0*pi/60.0; 
wr_ref=0 .2 * wr_rated; 

tmx=2.0*tl_rated; 

tmn=-2.0*tl_rated; 

clrscrQ; 

for(t=0 . 0 ;t<=5 0 . 0 ;t+=h) 

{ 

if(t>=25.0) 

wr_ref=-0 .2 * wr_rated ; 

if(t>=35.0) 

wr_reL=0.2*wr_rated; 
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/* SENSING OF ANALOG SIGNAL */ 


outportb(BASE+2,muxscan)'; 


for(i— 0;i<— 3 ;i++) 

{ 

outportb(BASE+0, trigger); 

while(l) 

{ 

u=inportb(BASE+8); 
check=u; 
check&= 0x80; 
if(! check) 
break; 

} 

Isb [i]=inportb(B ASE+0) ; 
msb[i]=inportb(BASE+l); 

} 

lsb[0]=lsb[0]»4; 

msb [0]=nisb [0] «4 ; 

val[0]=lsb[0]+msb[0]; 

val[0]=(val[0]-2047)*5/2047; 

Ia=val[0]*1.618; 

// printf("Ia=%2.if\t",Ia); 

lsb[l]=lsb[l]»4; 
msb[l]=msb[l]«4; 
val[l]=lsb[l]+msb[l]; 
val[l]=(val[l]-2047)*5/2047; 
Ib=val[l]* 1.678; 

// printf("Ib=%2.1f\n", Ib); 

lsb[2]=lsb[2]»4; 
msb[2] =msb [2] «4 ; 
val[2]=lsb[2]+msb[2]; 
val[2]=(val[2]-2047)*5/2047; 
vdc=val[2]* 148.22; 

// prmtf("vdc=%2.1f\n", vdc); 

lsb[3]=lsb[3]»4; 



insb[3]=msb[3]«4; 
val[3]=lsb[3]+msb[3]; 
val[3]=(val[3]-2047)* 5/2047; 
wr=val[3]*377.47*(2*pi/30.0)*0.441; 

// printf("wr=%2.1f\n", wr); 

/* END OF SENSING SIGNALS */ 

Ic=-Ia-Ib; 

ids=Ia; 

iqs=(Ib-Ic)/sqrt(3 .0); 
wre=wr_ref-wr; 

te_ref=te_ref+kp*(wre-wreO)+ki*h*wre; 

if(te_re£>=tmx) 

te_ref=tmx; 

if(te_ref<=tmn) 

te_ref=tmn; 

wre0=wre; 

ted=te_ref-te; 

if(ted>=0.0 && ted<htl && ct==0) 
ct=0; 

else if(ted>=htl && ted <ht2 && ct==0) 
ct=l; 

else if(ted>=ht2) 
ct=2; 

else if(ted>=htl && ted<ht2 && ct==2) 
ct=2; 

else if(ted<htl && ted>=0.0 &&ct=2) 
ct=l; 

else if(ted<=0.0 && ct==l) 
ct=0; 

else if(ted<=0.0 && ted>-htl && ct==0) 
ct=0; 

else if(ted <= -htl && ted>-ht2 && ct~0) 
ct=-l; 

else if(ted<=-ht2) 
ct=-2; 

else if(ted<= -htl && ted>-ht2 && ct==-2) 
ct=-2; 

else if(ted>-htl && ted<=0.0 &&ct==-2) 
ct=-l; 

else if(ted<0.0 && ted>=-htl && ct==-l) 
ct=-l; 



else if(ted>=0.0 && ct==-l) 
ct=0; 


fh=fref+hf; 

fl=fref-hf; 

if(fs<=fl) 

cf=l; 

if(fs>=fh) 

cf=0; 


// Calculation of flux vector 

if(tetae>=0.0 && tetae<pi/6.0) 
tetan=l; 

else if(tetae>=pi/6.0 && tetae<pi/3.0) 
tetan=2; 

else if(tetae>=pi/3.0 && tetae<pi/2.0) 
tetan=3; 

else if(tetae>=pi/2.0 && tetae<2.0*pi/3.0) 
tetan=4; 

else if(tetae>=2.0* pi/3.0 && tetae<5.0*pi/6.0) 
tetan=5; 

else if(tetae>=5.0*pi/6.0 && tetae<pi) 
tetan=6; 

else if(tetae>=pi && tetae<7.0* pi/6.0) 
tetan=7; 

else if(tetae>=7.0*pi/6.0 && tetae<4.0*pi/3.0) 
tetan=8; 

else if(tetae>=4.0*pi/3.0 && tetae<3.0*pi/2.0) 
tetan=9; 

else if(tetae>=3.0*pi/2.0 && tetae<10.0*pi/6.0) 
tetan=10; 

else if(tetae>=10.0*pi/6.0 && tetae<11.0*pi/6.0) 

tetan=ll; 

else tetan=12; 

// Selection of voltage vector from look up table 

v=lookuptable[cf] [ct+2] [tetan- 1 ] ; 
switch(v) 

{ 

case 1 1 : 

vds=(1.0/3.0)*vdc; 

vqs=0.0; 

sa 1 =0 ; sa2=0 ;sb 1 =0 ; sb2= 1 ;sc 1 =0 ;sc2= 1 ; 
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break; 
case 12: 

vds=(2.0/3.0)*vdc; 

vqs=0.0; 

sal =0;sa2=0;sb 1 = 1 ;sb2=0;sc 1 = 1 ;sc2=0; 
break; 
case 2 1 : 

vds=(1.0/6.0)*vdc; 
vqs=vdc/(2.0* 1.732); 
sal=0;sa2=0;sbl=0;sb2=0;scl=0;sc2=l; 
break; 
case 22: 

vds=(1.0/3.0)’'‘vdc; 

vqs=vdc/(1.732); 

sal =0;sa2=0;sb 1 =0;sb2=0;sc 1 = 1 ;sc2=0; 
break; 
case 3 1 : 

vds=-(1.0/6.0)*vdc; 
vqs=vdc/(2.0* 1.732); 
sa 1 =0;sa2= 1 ;sb 1 =();sb2=0;sc 1 =0;sc2= 1 ; 
break; 
case 32: 

vds=-(1.0/3.0)*vdc; 

vqs=vdc/1.732; 

sa 1 = 1 ; sa2=0 ;sb 1=0; sb2=0 ; sc 1 = 1 ;sc2=0 ; 
break; 
case 41: 

vds=-(1.0/3.0)*vdc; 

vqs=0.0; 

sal=0;sa2=l;sbl=0;sb2=0;scl=0;sc2=0; 
brealc; 
case 42: 

vds=-(2.0/3.0)*vdc; 

vqs=0.0; 

sal = 1 ; sa2=0 ; sb 1 =0 ; sb2=0 ; sc 1 =0 ;sc2=0 ; 
break; 
case 51: 

vds=-(l .0/6.0)*vdc; 
vqs=-vdc/(2.0* 1.732); 
sa 1 =0;sa2= 1 ;sb 1 =0;sb2= 1 ;sc 1 =0;sc2=0; 
break; 
case 52: 

vds=-(1.0/3.0)*vdc; 

vqs=-vdc/(1.732); 

sal=l;sa2=0;sbl=l;sb2=0;scl=0;sc2=0; 

break; 


case 61: 

vds=(l. 0/6.0)* vdc; 
vqs=-vdc/(2.0* 1.732); 
sal=0;sa2=0;sbl=0;sb2~l;scl=0;sc2=0; 
break; 
case 62: 

vds=(l .0/3.0)* vdc; 
vqs=-vdc/1.732; 

sal=0;sa2=0;sbl=l;sb2=0;scl=0;sc2=0; 
break; 
case 13: 

vds=(1.0/2.0)*vdc; 
vqs=vdc/(2.0* 1.732); 
sal =0;sa2=0;sb 1 =0;sb2= 1 ;sc 1=1 ;sc2=0; 
break; 
case 23: 

vds=0.0; 

vqs=vdc/(1.732); 

sal =0;sa2=l ;sb 1 =0;sb2=0;sc 1 = 1 ;sc2=0; 
break; 
case 33: 

vds=-(l .0/2.0)*vdc; 
vqs=vdc/(2.0* 1.732); 
sa 1 = 1 ; sa2=0 ; sb 1 =0 ;sb2=0 ;sc 1 =0 ; sc2= 1 ; 
break; 
case 43: 

vds=-(1.0/2.0)*vdc; 
vqs=-vdc/(2.0* 1.732); 
sal = 1 ;sa2=0;sb 1 =0;sb2= 1 ;sc 1 =0;sc2=0; 
break; 
case 53: 

vds=0.0; 

vqs=-vdc/(1.732); 

sal=0;sa2=l;sbl=l;sb2=0;scl=0;sc2=0; 
break; 
case 63 : 

vds=(l. 0/2.0)* vdc; 
vqs=-vdc/(2.0* 1 .732); 
sal=0;sa2=0;sbl=l;sb2=0;scl=0;sc2=l; 
break; 
case 7: 

vds=0.0; 

vqs=0.0; 

sal=l;sa2=l;sbl=l;sb2=l;scl=l;sc2=l; 

break; 


case 8: 
vds=0.0; 
vqs=0.0; 

sal=l ;sa2=l ;sbl=l ;sb2=l ;scl=l ;sc2=l ; 
break; 
default: 

vds=0.0; 

vqs=0.0; 

sal=l ;sa2=l ;sbl=l ;sb2=l ;scl=l ;sc2=l ; 
break; 

} 

fqs 1 =ls * iqs+lm* iqr; 
fds 1 =ls* ids+lm* idr; 
fqr 1 =lr * iqr+lm* iqs ; 
fdr 1 =lr*idr+Im*ids; 
wrl=wr; 


for(i=0;i<4;++i) 

{ 

kl [i]=h*(vqs-(rs/lls)*(fqsl-((lml/lls)*fqsl + (lml/llr)*fqrl))); 
k2 [i]=h* (vds-(rs/lls) * (fds 1 -(Qm 1 /11s)* fds 1 + (Im 1 /Ur) * fdr 1 ))) ; 
k3[i]=h*(wrl*fdrl-(rr/llr)*(fqrl-((lml/lls)*fqsl + (lml/llr)*fqrl))) ; 
k4[i]=h*(-wrl*fqrl-(rr/llr)*(fdrl-((lml/lls)*fdsl + (lml/llr)*fdrl))) ; 


if(i==0||i==l) 

{ 

fqsl+=kl[i]/2.0; 

fdsl+=k2[i]/2.0; 

fqrl+=k3[i]/2.0; 

fdrl+=k4[i]/2.0; 

} 

if(i==2) 

{ 

fqsl+=kl[i]; 

fdsl+=k2[i]; 

fqrl+=k3[i]; 

fdrl+=k4[i]; 

} 

} 

fqs+=(kl [0]+2.0*kl [l]+2.0*kl [2]+kl [3])/6.0; 
fds+=(k2[0]+2.0*k2[l]+2.0*k2[2]+k2[3])/6.0; 
fqr+=(k3[0]+2.0*k3[l]+2.0*k3[2]+k3[3])/6.0; 
fdr+=(k4[0]+2.0*k4[l]+2.0*k4[2]+k4[3])/6.0; 



fs=sqrt(fds * fds+fqs * fqs) ; 


fqm=((lml/lls)*fqs + (Iml/llr)*fqr); 
fdm=((lml/lls)*fds + (lml/llr)*fdr); 

iqr=(fqr-fqm)/llr; 

idr=(fdr-fdm)/llr; 

te=(3 .0/2.0)* (p/2 . 0) * (fds *iqs-fqs * ids) ; 
is=sqrt(ids*ids+iqs*iqs); 

if(fqs==0.0&&fds==0.0) 

tetae=0.0; 

else 

tetae=atan2(fqs,fds); 
if(tetae<0.0) 
tetae= tetae+2.0*pi; 

/^.KHc********** SWITCHING OUT CONTROL puLSES*************/ 

swlogic = 32 * scl + 16 * sc2 + 8 * sbl + 4 * sb2 + 2 * sal + sa2; 
outportb(BASE+3,swlogic); 

/** SENDING ANALOG SIGNALS OUT*************/ 
vad=2047.0+2047.0 *fqs/25.; 
vad=floor(vad); 
mb=vad/16; 
lb=vad-inb*16; 
lb=lb«4; 

outportb(BASE+6,lb); 

outportb(BASE+7,mb); 

vad=2047.0+2047. * te/15.0; 

vad=floor(vad); 

mb=vad/16; 

lb=vad-mb*16; 

lb=lb«4; 

outportb(BASE+4,lb); 
outportb(B ASE+5 ,mb) ; 

} 

end=clock(); 

printf("%f\n",((end-strt)/CLK_TCK)*h); 

// fclose(fpl); 

} 
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